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ABSTRACT 



■ Context. The XMM-Newton Distant Cluster Project is a serendipitous survey for clusters of galaxies at redshifts z > 0.8 based on deep archival 

■ XMM-Newton observations. X-ray sources identified as extended are screened against existing optical all-sky surveys for galaxies, in case of 
i | ' candidate high-z clusters followed up with imaging at 4m-class telescopes and, ultimately, multi-object spectroscopy at 8m-class telescopes. Low- 

, significance candidate high-z clusters are followed up with the seven-channel imager GROND (Gamma-Ray Burst Optical and Near-Infrared 
Detector) that is mounted at a 2m-class telescope. Its unique capability of simultaneous imaging in the g', r', i', z', J, H, Ks bands enables the use 
of the photometric redshift technique. 

pH ' Aims. Observing strategy, data reduction and analysis, depth and accuracy of the simultaneous multi-wavelength photometry are discussed with 

the goal of establishing GROND as a useful instrument to confirm X-ray selected (high-z) clusters. 
I . Methods. The test case is XMMU J0338. 7+0030, suggested to be at z ~ 1.45 ± 0.15 (lcr) from the analysis of the z - H vs H colour-magnitude 
O ' diagram obtained from the follow-up imaging. Later VLT-FORS2 spectroscopy enabled us to identify four members, which set this cluster at 
H ' Z = 1-097 ± 0.002 (lcr). To reach a better knowledge of its galaxy population, we observed XMMU J0338.7+0030 with GROND for about 6 hr. 
C/3 The publicly available photo-z code le Phare was used. 

■ Results. The Ks-band number counts of the non-stellar sources out of the 832 detected down to z' ~ 26 AB mag (lcr) in the 3.9 x 4.3 arcmin 2 
1 1 ■ region of XMMU J0338. 7+0030 imaged at all GROND bands clearly exceed those computed in deep fields/survey areas at ~ 20.5-22.5 AB mag. 

The photo-z's of the three imaged spectroscopic members yield z = 1.12 ± 0.09 (lcr). The spatial distribution and the properties of the GROND 
sources with a photo-z in the range 1.01-1.23 confirm the correspondence of the X-ray source with a galaxy over-density at a significance of at least 
4.3 cr. Candidate members that are spectro-photometrically classified as elliptical galaxies define a red locus in the V — z' vs z' colour-magnitude 
diagram that is consistent with the red sequence of the cluster RDCS J0910+5422 at z = 1.106. XMMU J0338.7+0030 hosts also a population of 
bluer late-type spirals and irregulars. The starbursts among the photometric members populate both loci, consistently with previous results. 
Conclusions. The analysis of the available data set indicates that XMMU J0338. 7+0030 is a low-mass cluster (M 2 ao ~ 10 14 M©) at z = 1.1. 
With the photometric accuracy yielded by the present unoptimized multi-band observations with GROND, we not only confirm the spectroscopic 
redshift of this cluster but also show that it hosts a galaxy population that can still undergo significant bursts of star-formation activity. 
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. ^ | 1 . Introduction hot plasma^ emits), or at sub-millimeter wavelengths (owing to 

. the Sunyaev-Zel'dovich - SZ - effect, Sunyaev & Zel'dovich 

5^ ■ A process of self-similar gravitational clustering in an expand- 1972). In particular, the X-ray or SZ-effect detection ensures that 

5^ \ ing universe leads to the build-up of structures in a hierarchical a galaxy system is truly bound and associated with a deep grav- 

fashion across time (White & Rees 1978; White et al. 1987). itational potential well. X-ray observations of groups/clusters 

Ensembles of galaxies and hot (T > 10 6 K) plasma within gravi- can now be pushed towards the distant universe (i.e., z > 0.8) 

tationally bound haloes dominated by dark matter define groups thanks to the Chandra and XMM-Newton space-borne observa- 

and clusters, with total masses from ~ 10 13 to ~ 10 15 M©, re- tories (e.g., Stanford et al. 2001, 2002; Rosati et al. 2004; Mullis 

spectively. These objects offer a biased view of the evolution of e t al. 2005; Finoguenov et al. 2007; Suhada et al. 2010; however, 

the large-scale structure (Kaiser 1984) that is easy to model the- see Rosati et al. 1998). With the SZ-effect this has become pos- 

oretically from ab-initio principles. sible very recently but only for massive clusters (Staniszewski et 

The identification of a group/cluster requires observations al - 2009 '- Vanderlinde et al. 2010; Marriage et al. 201 1). 

at rest-frame optical/near-infrared (IR) wavelengths (where the Compiling a sample of distant galaxy clusters has a twofold 

photospheric stellar emission dominates), in X-rays (where the importance. Firstly, it offers a way to trace structure growth in 

the far universe that enables the study of the effects of cold dark 



Whether the hot X-ray emitting plasma belongs to a group or a 



Currently visiting astronomer at the MPE cluster, we will refer to it as the intracluster medium (ICM). 
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matter (CDM) and dark energy (DE) on the cosmic structure 
evolution (e.g., Albrecht et al. 2006; Vikhlinin et al. 2009; Mantz 
et al. 2010a). Secondly, it fosters the investigation of the evo- 
lution of the galaxy population, the baryon mass component of 
groups/clusters, the thermodynamics and chemical abundance of 
the ICM (e.g., Postman et al. 2005; Mei et al. 2006a,b; Balestra 
et al. 2007; Gobat et al. 2008; Maughan et al. 2008; Giodini et 
al. 2009; Mantz et al. 2010b; Santos et al. 2010; Strazzullo et al. 
2010). About twenty clusters at z = 0.8-2.1 have been discov- 
ered and studied in some detail (e.g., Stanford et al. 2001, 2002, 
2006; Blakeslee et al. 2003; Mullis et al. 2005; Demarco et al. 
2007; Hilton et al. 2007, 2009; Eisenhardt et al. 2008; Lamer et 
al. 2008; Menci et al. 2008; Wilson et al. 2009; Kurk et al. 2009; 
Tanaka, Finoguenov & Ueda 2010; Papovich et al. 2010; Henry 
etal. 2010; Gobat et al. 2011;Fassbenderet al. 201 la). However, 
a better sampling in redshift and total mass is necessary to build 
a representative sample of high-z groups and clusters. Existing 
X-ray surveys allow clusters in the redshift range 0.8-1.5 to be 
detected, with about 1 .5 such clusters per square degree down to 
a flux level of ~ 6 x 10~ 15 erg s" 1 cm" 2 (Suhada et al. 2012). 

The XMM-Newton Distant Cluster Project (XDCP) is a 
serendipitous survey for galaxy clusters at z > 0.8 based on deep 
archival XMM-Newton observations: it now covers an area of 
80 deg 2 (Bohringer et al. 2005; Fassbender 2007; Fassbender et 
al. 201 lc). XDCP is concerned with the build-up of a represen- 
tative and complete sample of distant clusters and the accurate 
determination of the total mass of each system. Both aspects are 
crucial for a test of the current understanding of structure for- 
mation and evolution in ACDM universes (e.g., Hoyle, Jimenez 
& Verde 2011). In XDCP, X-ray sources identified as extended 
are screened against existing optical all-sky surveys and even- 
tually followed up with imaging in the z, H or I, J bands at 
4m-class telescopes (occasionally in the r, Z bands at 8m-class 
telescopes). If the distribution of the galaxies associated with 
an X-ray source exhibits evidence of a so-called "red sequence" 
(e.g. Bower, Lucey & Ellis 1992; Kodama et al. 1998; Gladders 
& Yee 2000; see Arimoto & Yoshii 1987 for its physical in- 
terpretation) in the ensuing colour-magnitude diagram, a can- 
didate high-z cluster is robustly selected and proposed for spec- 
troscopic confirmation at 8m-class telescopes. This strategy has 
already yielded more than 30 spectroscopically confirmed dis- 
tant clusters since the discovery of XMMU J2235.3-2557 at z — 
1.393 (Mullis et al. 2005). They include XMMU J03 3 8. 8 +0021 
at z = 1.49 (Nastasi et al. 2011), XMMU J1007.4+1237 at 
z = 1.555 (Fassbender et al. 2011b) and XMMU J0044.0-2033 
atz= 1.579 (Santos etal. 2011). 

Among the main results from XDCP we quote the multi- 
wavelength studies of XMMU J1229+0151 atz = 0.975 (Santos 
et al. 2009) and XMMU J2235.3-2557 (Rosati et al. 2009; 
Strazzullo et al. 2010); the weak-lensing estimate of the total 
mass of XMMUJ2235.3-2557 (Jee et al. 2009); the total mass 
estimate of XMMU J100750.5+125818 (z = 1.082) from strong 
lensing, optical spectroscopy and X-ray imaging (Schwope et 
al. 2010); the pan-chromatic study of XMMU J1230.3+1339 at 
z = 0.975 (Fassbender et al. 2011a) and its total mass estimate 
from weak-lensing analysis of ground-based imaging (Lerchster 
et al. 2011); the evidence of starburst activity in the core of 
XMMU J1007.4+1237 (Fassbender et al. 201 lb). 

In order to maximize the numerical throughput from XDCP, 
low-significance candidate high-z clusters are followed up 
with the seven-channel imager GROND (Gamma-Ray Burst 
Optical and Near-Infrared Detector, Greiner et al. 2008), that 
is mounted at a 2m-class telescope, before being proposed for 
multi-object spectroscopy (MOS) at 8m-class telescopes. The 



optical confirmation of these distant clusters lies in the discov- 
ery of a significant feature in the distribution of the photometric 
redshifts of galaxies selected in a region corresponding to the 
extended X-ray emission. Otherwise, it relies upon the fact that 
colour- or photo-z-selected galaxies within a given projected dis- 
tance from the centroid of the extended X-ray emission define an 
over-density with respect to the surrounding region and/or a red 
sequence in a given colour-magnitude diagram. 

Here we report on the feasibility of this approach by il- 
lustrating results obtained from the multi-wavelength obser- 
vations of XMMUJ0338.7+0030. Follow-up z,H photometry 
was originally obtained with the near-IR, wide-field camera 
OMEGA2000 mounted at the prime focus of the 3.5m tele- 
scope at Calar Alto, Spain. A redshift z ~ 1.45 + 0.15 was 
inferred from comparison between the average z - H colour 
of the reddest galaxies within 45" from the X-ray position 
of XMMUJ0338.7+0030 and the synthetic z - H colours ex- 
pected for a large suite of simple stellar population (SSP) mod- 
els (see Fassbender et al. 201 lc). The synthetic photometry was 
computed using the spectro-photometric evolutionary synthesis 
model PEGASE v2 (Fioc & Rocca-Volmerange 1997) in analogy 
with Pierini et al. (2005) and Wilman et al. (2008). A calibration 
of the best-performing SSP model was originally based on ob- 
served red-sequence galaxies in 10 X-ray selected clusters with 
spectroscopic redshifts up to 1.5 (Fassbender 2007). This fidu- 
cial model corresponds to a Salpeter (1955) stellar initial mass 
function (IMF) between 0.1 and 120 Mq (as all the others), a 
solar metallicity and a formation redshift equal to 5. 

In addition to its being potentially one of the few clusters at 
z > 1.4 known so far, XMMU J0338.7+0030 represented a test 
case for XDCP sources at the detection limit of the survey. It was 
an ideal target for GROND imaging because its region is covered 
by the Sloan Digital Sky Survey (SDSS; York et al. 2000) and the 
2 Micron All Sky Survey (2MASS; Skrutskie et al. 2006), which 
enables a robust photometric calibration for all seven bands. For 
the same region, sparse redshift information was later available 
from observations with the spectrograph FORS2 mounted at the 
ESO Very Large Telescope (VLT) at Paranal, Chile. 

The XMM-Newton, OMEGA2000, GROND and FORS2 
observations of XMMU J0338.7+0030, plus the reduction and 
analysis of the ensuing data, are discussed in Sect. 2. Estimates 
of the X-ray bolometric luminosity, temperature and total mass 
of this cluster, plus the spectroscopic determination of its red- 
shift, are given there. Section 3 contains results obtained from 
the GROND data, i.e. galaxy number counts, the photo-z distri- 
bution of the detected sources, the two-dimensional (2-D) spatial 
distribution of photo-z-selected galaxies and their distribution 
in specific colour-magnitude diagrams. Discussion and conclu- 
sions appear in Sect. 4 and 5, respectively. 

Throughout we adopt a ACDM cosmological model (O m = 
0.3, Q A = 0.7) with H = 70 kms~' Mpc -1 , that is broadly 
consistent with the main results from the 5 yr operations of the 
Wilkinson Microwave Anisotropy Probe (WMAP5, Dunkley et 
al. 2009; Komatsu et al. 2009). Accordingly, at the spectroscopic 
redshift of XMMU J0338.7+0030 (z = 1.097; Sect. 2.4.2), the 
age of the Universe is 5.4 Gyr, the cosmic evolution factor 
E(z) = H(z)/H Q is equal to 1.86, the luminosity distance is 
7409.2 Mpc, the angular scale is 8.168 kpc/". Hence, an angular 
distance of 1' corresponds to ~ 490 kpc at the cluster distance. 
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Fig. 1. A z+H colour composite image (4.5' on the side) of the environment of XMMU J0338. 7+0030 with an overlaid contour map 
of the X-ray emission detected by XMM-Afewfon (in yellow). These log-spaced X-ray contours correspond to significance levels of 
2-16cr; they were derived from the adaptively smoothed, combined images of the cluster environment. Red contours highlight the 
over-density of galaxies with very red z - H colours (3.0 < z - H < 4.8 Vega mag), as obtained from the OMEGA2000 imaging (see 
Sect. 2.2.2). White dashed and solid lines encircle regions within 30" and 60" from the X-ray centroid of the cluster (R. A.(J2000.0): 
03h38m44.2s, Dec(J2000.0): +00d30m01.8s), respectively; they guide the eye. The red contour closest to the X-ray position of the 
cluster corresponds to a significance of 4.6<x above the mean background. 



2. Observations, data reduction and analysis 

2.1. X-ray imaging with XMM-Newton 
2.1.1. Initial X-ray selection 

A total of 470 XMM-Newton (Jansen et al. 2001) archival 
fields with nominal exposure times longer than 10 ks have 
so far been processed and analyzed in XDCP as exten- 
sively described in Fassbender (2007) and Fassbender et al. 
(2011c). The serendipitous, weak X-ray source identified as 
XMMUJ0338.7+0030 was detected in the field with observa- 
tion identification number (OBSID) 0036540101. The target 
was SDSS 033829.31+002156.3 (R.A.(J2000.0): 03h38m38.5s, 
Dec(J2000.0): +00d39m28s), a quasar at z = 5.07091 (Peroux et 
al. 2001) imaged on February 22, 2002 with XMM-Newton for 
a nominal exposure time of 22.9 ks. 

This serendipitous X-ray source was found at an off- 
axis angle of 9.6' (R.A.(J2000.0): 03h38m44.2s, Dec(J2000.0): 
+00d30m01.8s) and a significance of about 5<x during the early 
XDCP source detection run performed with SAS v6. 1 . The tasks 
eboxdetect was used for a sliding box detection, followed by a 
subsequent maximum likelihood analysis with emldetect. This 
yielded an extent of 3.48 pixels (i.e., 13.92 arcsec), a source ex- 
tent likelyhood of about 2<x and a flux of ~ 8 x 1 0~ 15 erg s~ 1 cm -2 
for ~ 32 PN counts, close to the detection threshold. Hence, 
XMMU J0338. 7+0030 appeared as a marginally extended X-ray 
source at the XMM-Newton detection threshold for extended 
sources. An X-ray contour map is overlaid on the z+H colour 
composite image obtained from the medium-deep imaging with 
OMEGA2000 (Sect. 2.2.1) in Fig. 1. 



The circular area centred on the X-ray emission centroid of 
XMMUJ0338.7+0030 and delimited by a radius of 1' is as- 
sumed as the "bona fide" cluster region throughout. 

2.1.2. The upgraded X-ray picture of XMMU J0338.7+0030 

The XMM-Newton field containing XMMU J0338.7+0030 was 
reprocessed with SASvlO.0.0 following Suhada et al. (2011). 
These observations are strongly flared and an automatic flare re- 
moval is impossible; hence, the contaminated periods were ex- 
cised manually and the automatic two-step flare cleaning process 
was run on the remaining good part of the data for the removal 
of high background periods. The loss of information is severe: 
clean net exposure times on axis amount to 9 ks for either MOS 
camera and 7.8 ks for the PN instrument, whereas the effective 
(vignetted) clean net exposure times on source amount to 4.5- 
5.5 ks for all detectors. Even after cleaning a slight residual qui- 
escent contamination (about 10% above normal) remains; large 
part of it should be captured by the background model, so this 
anomalous contamination likely has a negligible effect. 

With SASvlO.0.0 XMMU J0338.7+0030 is detected as a 
very weak source without a significant extent, located about 10" 
to the east with respect to the X-ray position originally estimated 
(Sect. 2. 1 . 1). A Kolmogorov-Smirnov (KS) test with the null hy- 
pothesis that the extracted X-ray profile corresponds to that of a 
point-like source gives -log Pks = 2.08 and -log Pks = 1-57 
for the PN and MOS data, respectively. In other words, there is 
a chance of 0.8% or 2.7% that XMMU J0338.7+0030 is a point- 
like source. The result based on the PN instrument is likely un- 
realistic, because PN has a chip gap very close to the centre. 
Furthermore, the performed test uses the on-axis point-spread- 
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function (PSF) of XMM-Newton, so the output probabilities 
are slightly underestimated for the off-axis source identified as 
XMMU J0338. 7+0030. We conclude that this source is weak 
and marginally extended. 

This is supported by the regular behaviour of the growth 
curves, where the growth curve analysis (GCA) method of 
Bohringer et al. (2000) was applied to measure the flux in the 
soft 0.5-2 keV band. The two growth curves obtained with the 
MOS cameras or the PN instrument agree (see Fig. 2). However, 
they are not independent: the latter was slightly corrected to 
obtain a rough agreement, because PN has a chip gap very 
close to the centre as already mentioned. The double compo- 
nent background model was used in the first instance; the alter- 
native use of a spline fitting led to the same results. As a re- 
sult, we measured a flux /x,500 — 7.1(±2.3) x 10~ 15 erg s _1 cm" 2 
within a circular aperture of radius equal to /?5oJ3 (correspond- 
ing to 51.2"). At the cluster redshift, this translates into an X- 
ray luminosity of Lx.soo — 4.5(±1.4) x 10 43 erg s" 1 in the soft 
band or L^ 0Q =* l'.l(±0.3) x 10 44 ergs" 1 for the bolomet- 
ric luminosity. The latter quantity was inferred as described in 
Suhada et al. (2010). In brief, the growth curve was used to iter- 
atively obtain a self-consistent set of parameters using the X-ray 
luminosity-temperature and X-ray luminosity-total mass rela- 
tions from Pratt et al. (2009). 

Owing to the faintness of the source with less than 50 net 
counts in the soft band, the determination of additional struc- 
tural or spectral parameters is currently not feasible. The po- 
tential flux contribution of point sources can presently not be 
accounted for either, which implies that the stated luminosities 
have to be taken as upper limits. However, based on the lumi- 
nosity scaling relation in Pratt et al. (2009), first rough estimates 
on the expected temperature of Tx ~ 2.4 keV and the total mass 
of M 500 ~ 7 x 10 13 M (i.e., M 20 o ~ 10 14 M ) can be obtained 
from the previous values. Hence, XMMU J0338. 7+0030 appears 
to be a system at the boundary between X-ray selected groups 
and clusters, as XMMU J1532.2-0836 at z = 1.358 (Suhada et 
al. 2011) and XMMU J0338.8+0021 at z = 1.49 (Nastasi et al. 
2011). 

An exposure-corrected count image of the innermost region 
of the cluster is produced after masking point-like sources (cf. 
Fig. 1) and is then adaptively smoothed with a circular top hat 
filter. The ensuing X-ray contours are reproduced in the right 
panel of Fig. 4 (Sect. 2.3.1) using a squared scale. 

2.2. Near-infrared follow-up imaging with OMEGA2000 
2.2.1. Observations and data reduction 

XMMUJ0338.7+0030 was observed with the near-IR camera 
OMEGA2000 (Bailer-Jones, Bizenberger & Storz 2000) that 
is mounted at the 3.5m telescope at the Centra Astronomico 
Hispano Aleman (CAHA) at Calar Alto, Spain. OMEGA2000 
provides a field-of-view (FoV) of 15.4' x 15.4'. In particular, 
medium-deep imaging of the XMMU J0338. 7+0030 region in 
the H-band (50 min) and z-band (53 min) was obtained on 
January 5, 2006 under good but non-photometric conditions 
of the sky and with a seeing (or PSF) that varied between 
1.0" and 1.5" (z-band) or 1.6" (H-band) in full-width-at-half- 
maximum (FWHM). The same region was reobserved in z-band 
(5 min) on October 30, 2006 in photometric conditions of the 



2 Ra (A = 500) is the radius within which the total mass density of a 
group/cluster is equal to A times the critical density of the universe (p c ). 
Correspondingly, Ma = A p c (z) (4^/3) R 3 A is the total mass within ^a- 
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Fig. 2. Results from the application of the GCA method of 
Bohringer et al. (2000) in the soft 0.5-2 keV band to the MOS 
and PN observations of XMMU J0338.7+0030. The two growth 
curves are well behaved and agree with each other. 

sky. Observations of a designated SDSS standard star (Smith et 
al. 2002) were performed then to enable a photometric calibra- 
tion of the science frames. 

These data were reduced with the OMEGA2000 near-IR 
pipeline (Fassbender 2007), which performs the basic frame cal- 
ibrations (flat-fielding, bad pixel correction), a two-pass back- 
ground subtraction with masked objects, frame alignment and 
optimally weighted co-addition of the individual 40 s (60 s) ex- 
posures in H-band (z-band). For each band, the individual frames 
were stacked without PSF-matching, which gave a measured 
seeing of 1.3" (z-band) or 1.4" (H-band) in the final image. 

For the CAHA data, source detection, extraction of pho- 
tometry and star/galaxy separation are described in detail in 
Fassbdender et al. (2007). The SExtractor (Bertin & Arnouts 
1996) photometry was calibrated to the Vega magnitude system 
using 2MASS point sources (Cutri et al. 2003) in H-band and 
designated SDSS standard star observations in z-band, cross- 
checked with SDSS photometry in the science frame. In partic- 
ular, total magnitudes were measured in a standard Kron ( 1980) 
aperture of factor 2.5 and minimum radius of 3.5. The limit- 
ing (Vega) magnitude^] were determined to be H\i m ~ 21.2 and 
Zum ~ 23.1. They correspond to the expected apparent magni- 
tudes of an old, passively evolving galaxy at the estimated red- 
shift of z ~ 1-45 (see Sect. 1) for an absolute magnitude equal 
to ~ M* + 1.3 in H-band and ~ M* in z-band, where M* is the 
characteristic magnitude of the corresponding galaxy luminosity 
function at z = 1.45 (see Fassbender 2007). 

2.2.2. An over-density of very red galaxy at the location of 
XMMUJ0338.7+0030 

In Fig. 1 (Sect. 2.1.1), the projected distribution of galaxies with 
very red colours (3.0 < z - H < 4.8 Vega mag) exhibits a signif- 
icant over-density: 4.6<x above the mean background. The peak 
of this over-density is offest to the East by about 20" with respect 
to the X-ray position of XMMU J0338.7+0030. When consider- 
ing the degradation of the angular resolution of XMM-Newton 
(~ 14" half-energy-width - HEW) at the off-axis angle of this 

3 These limiting magnitudes correspond to a 50% completeness, as 
estimated from the number counts obtained from the OMEGA2000 
imaging only, as in Fassbender et al. (201 la). 



4 



D. Pierini et al.: I. The galaxy population of XMMU J0338.7+0030 at z = 1.1 



X-ray source (Sect. 2.1.1), an association between the previous 
over-density of red galaxies and the cluster is very likely. 

Figure 3 shows the z - H vs H colour-magnitude diagram 
for the galaxies detected in the XMMU J0338.7+0030 region 
imaged with OMEGA2000. In particular, seven galaxies within 
45" from the X-ray position of the cluster exhibit 2.8 <z — H< 
3.8 Vega mag). A redshift z ~ 1.45 + 0.15 was initially inferred 
from comparison between the average z - H colour of these 
galaxies and the synthetic z-H colours expected for the fiducial 
SSP model described in Sect. 1. The assumption was that such a 
template described all those galaxies. 

2.3. Optical/near-infrared imaging with GROND 

2.3.1 . Simultaneous images of XMMU J0338. 7+0030 

A region of the sky centred on the XMM-Newton X-ray position 
of XMMUJ0338.7+0030 was imaged on four different nights 
(September 24, 2008, and January 27, 28, 29, 2009), in con- 
ditions of clear sky but variable seeing, with GROND. This is 
a seven-channel optical/near-IR imager primarily designed for 
rapid, simultaneous multi-wavelength observations of gamma- 
ray burst afterglows and built by the high-energy group of MPE 
in collaboration with the Landessternwarte Tautenburg and ESO. 
GROND was commissioned at the MPI/ESO 2.2m telescope at 
La Silla, Chile, in April 2007 and is operated as a Pi-instrument 
(for its comprehensive description see Greiner et al. 2008). 

For the optical channels, the incoming telescope beam is split 
using four dichroics designed in such a way that their transmis- 
sion functions are identical to those of the Sloan g', r', i', z' filter 
system, with the exception of the i'-band. Each of the four iden- 
tical optical CCDs has a FoV of 5.4 x 5.4 arcmin 2 and a scale of 
0. 158"/pixel. The near-IR part of GROND is designed as a focal 
reducer system yielding a FoV of 10 x 10 arcmin 2 with a scale 
of 0.6"/pixel. The GROND J, H, Ks dichroic filters are close to 
other near-IR filter systems. Photometric system transformations 
between the GROND g',r',i',z' bands and the Sloan filter sys- 
tem and between the GROND J, H, K s band s and the 2MASS 
filter system are given in Greiner et al. (20Q8J. 

A total of 10 observations in the seven channels was executed 
through the GROND standard observing block (OB) named 
"20min4TD". By design, this OB corresponds to four single 
369 s-long exposures in the g',r' channels and four 369 s-long 
exposures out of four pairs of sub-exposures in the i',z' chan- 
nels, each long exposure being taken at one out of four different 
telescope dither positions. In parallel, the three near-IR chan- 
nels are operated with 10 s integrations. At each telescope dither 
position, a six-position dither pattern is executed through a flip 
mirror (internal dithering) only for the Ks channel; five 10 s- 
long exposures per mirror position are executed for the J, H, Ks 
channels. This produces a total of 120 single 10 s-long exposures 
per OB for each near-IR channel. Given the impact of overheads 
(about 50% of the total exposure time of 4 hr for each J, H, Ks 
channel), the total observing time amounted to about 6 hr. 

Pixel- and gain-corrected, astrometrized, stacked images for 
each channel and per OB were produced using the data re- 
duction and photometry tools of the GROND pipeline (Kiipcii 
Yoldas et al. [20081 Kruhle r et a l. 120081 , based on standard 
IRA^PyRAF tools (Tody [19931 . Astrometry was achieved 
through the use of stars in common with the SDSS catalogue 

4 IRAF is distributed by the National Optical Astronomy Observa- 
tories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



(Abazajian et al. 120091 ) f or the optical bands and the 2MASS 
catalogue (Skrutskie et al. 12006 ) for the near-IR ones. For each 
image the FWHM of the average point-spread-function (PSF) 
was measured from the surface brightness profiles of randomly 
distributed non-saturated stars. Images corresponding to values 
of the average seeing FWHM higher than 1.5" were not consid- 
ered in the following analysis. For each band, all selected im- 
ages were convolved to the worst PSF (~ 1 .5" FWHM) and then 
stacked with a median algorithm. The total exposure times of the 
stacked g',r', i', z', J, H, Ks images are equal to 1 1808 s, 14760 s, 
14760 s, 14760 s, 12000 s, 10800 s, and 9600 s, respectively. The 
stacked near-IR images were regridded onto the same scale as 
the optical ones (0.158"/pixel). All seven stacked images were 
re-mapped to the same aspect with the z'-image as reference. 

Figure 4 reproduces a J,z',g' colour-composite (RGB) 
image of the region of XMMU J0338.7+0030 mapped with 
GROND (left panel) together with the deep z'-image (right 
panel). 

2.3.2. The photometric catalogue of XMMU J0338.7+0030: 
source extraction and photometry 

A separate stack of seven individual images in the z' channel, 
PSF-matched to the same angular resolution of 1.1" (FWHM), 
was used as a detection image with a total exposure time equal 
to 10332 s. This choice was the best compromise between image 
quality, photometric depth in the red/near-IR channels, robust- 
ness of detection and penalty in the source detection and extrac- 
tion of photometry associated with deblending. It also enabled a 
better, independent separation between extended and point-like 
sources along the accessible flux domain. 

Automatic source detection and extraction of photometry 
within an area of 3.9 x 4.3 arcmin 2 was accomplished through 
the publicly available, standard software SExtractor (Bertin & 
Arnouts 1996). A source was defined by the excess over a detec- 
tion threshold equal to the standard deviation (l<x) of the global 
image background in at least two contiguous pixels; in addition, 
a Gaussian filter with a FWHM of 1.1" was applied for detec- 
tion. Photometry of individual sources was extracted in six cir- 
cular apertures with diameters equal to 1,2, 3, 4, 5 and 10 arcsec, 
in a Rron (1980) aperture of factor 2.5 and minimum radius of 
3.5 (SExtractor MAG_AUTO) and in a Petrosian (119761 ) aperture 
of factor 2.0 and minimum radius of 3.5. 

In the first step of extraction of source photometry, typical 
photometric zeropoints for the GROND channel^ were assumed 
together with typical corrections for atmospheric absorption (or 
airmass), in order to aid identification of stars in common with 
the existing SDSS and 2MASS catalogues and enable the cali- 
bration of the zeropoints of the present g',r',i',z', J, H, Ks im- 
ages. In particular, we made use of the seventh Sloan data re- 
lease (DR7; Abazajian et al. 2009 and references therein) and 
the 2MASS All Sky Catalog of point sources (Cutri et al. 2003). 
After this photometric calibration, adjustments for the newly de- 
termined zeropoints were introduced and magnitudes were cor- 
rected for Galactic extinction according to the values tabulated 
in Schlegel et al. ( 1998 ) in addition to the previous correction for 
atmospheric absorption. 



5 GROND photometric zeropoints in the AB magnitude system are 
available at http://www.mpe.mpg.de/~jcg/GROND/calib.html and in 
Greiner et al. {2008 1 together with conversions between different filter 
systems or from the AB to the Vega magnitude system. These conver- 
sions are equal to -0.008, -0.151, -0.386, -0.515, -0.910, -1.381 and 
-1.795 mag for, respectively, the GROND g',r',i',z', J, H, Ks bands. 
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Fig. 3. z - H vs H colour-magnitude diagram for the region encompassing XMMU J0338. 7+0030 imaged with OMEGA2000. 
Red filled circles represent galaxies within 45" from the X-ray position of the cluster. The blue dashed lines reproduce the 50% 
completeness of the photometric catalogue extracted from the OMEGA2000 data. The horizontal red short-dashed line represents 
the z - H colour expected for an SSP model corresponding to the evolutionary phase of an instantaneous burst of star formation 
observed at z — 5 with Salpeter (1955) stellar IMF and solar metallicity observed at z — 1 .45. This SSP model describes the average 
colour of the seven galaxies associated with XMMU J0338. 7+0030 that exhibit 2.8 < z - H < 3.8 Vega mag. The vertical red 
dotted line represents the characteristic magnitude H* at z — 1 .45 (from Fassbender et al. 2007). Empty black squares mark the four 
spectroscopic cluster members that were successively established to be at z = 1.1 through VLT-FORS2 observations (see Fig. 1 1 
and Table 1 in Sect. 2.4.2); blue filled circles highlight those with 1D=12 and ID=15, which are at a cluster-centric distance greater 
than 45". The two spectroscopic cluster members without convincing evidence of ongoing star-formation activity (with ID=2 and 
ID=16) exhibit about as red z - H colours as the previous SSP model (see Appendix B, however). 



For the photometric calibration of the g',r',i',z' (J, H, Ks) 
images, the celestial coordinates of GROND sources and SDSS 
(2MASS) catalogue stars were first cross-correlated using a 
search radius of 3". "Machine" Petrosian magnitudes (corrected 
only for atmospheric absorption) and catalogue Petrosian mag- 
nitudes (not corrected for Galactic extinction) of non-saturated 
sources that are unambiguously identified as stars and have ro- 
bust photometry were contrasted after taking into account even- 
tual colour terms (see footnote 5). For a given band, the weighted 
mean of the magnitude differences was taken as the correction 
to the assumed zeropoint and the value of the weighted standard 
deviation as the uncertainty of the photometric zeropoint. These 
uncertainties are equal to ±0.009, +0.010, ±0.011, ±0.011, 
±0.041, ±0.009 and ±0.036 mag in the g', r', i', z', J, H, Ks chan- 
nels of GROND, respectively. They reflect the quality of the ref- 
erence photometry and the poor statistics of the matched stars 
used as calibrators. The good agreement between our set of cali- 
brated magnitudes and that of SDSS (2MASS) for all the 38 (9) 
stars in common is shown in Fig. 5 (6). 

Our calibrated Petrosian magnitudes are compared to the 
SDSS ones (all not corrected for Galactic extinction) for the 65 
sources in our photometric catalogue that are associated with 
galaxies identified in DR7 using a search radius of 3". In addi- 
tion, twenty-four of these 65 SDSS galaxies have a photometric 
redshift in DR7. The overall agreement between the two sets 
of Petrosian magnitudes is satisfying (Fig. 5). Unfortunately, no 
galaxy counterparts to our sources seem to exist in any 2MASS 
catalogue, according to NED. This is no surprise given the shal- 



low limiting magnitudes of 2MASS (see Jarrett et al. 2000; 
Skrutskie et al. 2006). 

A calibrated photometric catalogue with 832 sources was 
built as previously described, where magnitudes are corrected 
for atmospheric extinction and Galactic extinction. Limiting 
magnitudes (lcr values) were computed as the median values 
of magnitudes obtained in 3" apertures on 10 separated regions 
of blank sky distributed across the individual g', r', i', z', J, H, Ks 
images. They correspond to 27.30, 26.89, 26.54, 26.06, 24.40, 
23.50 and 22.91 mag, respectively. We note that the correspond- 
ing 5<x values in the z',J,H, Ks channels are about the depths 
in analogous bands that the ESO public survey VIKING ( VISTA 
Kilo-Degree Infrared Galaxy Survey; P.I.: W. Sutherland) is ex- 
pected to reach (see e.g., Arnaboldi et al. 2007). 

Corrections to the SExtractor output photometric errors for 
a given band (in part caused by correlated noise in the stacked 
images) were computed from the ratio of the median value of 
the uncertainties of the previous 10 aperture magnitudes and the 
median SExtractor uncertainty of the 3" aperture photometry of 
sources with magnitudes within ±lcr from the limiting magni- 
tude. These correction factors amount to 1.74, 1.67, 1.82, 1.28, 
2.15, 1.36, and 1.78 for the g',r',i',z',J, H, Ks bands, respec- 
tively. 

Extended and point-like sources were initially identified on 
the basis of surface brightness distribution (i.e., the SExtractor 
parameters stellarity index and flux radius) and apparent bright- 
ness. The photometry of point-like sources enabled us to esti- 
mate corrections from the 3" aperture magnitudes to magnitudes 
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Fig. 4. J, z', g' colour-composite image of the 3.9 x 4.3 arcmin 2 region of XMMU J0338. 7+0030 mapped with GROND at all seven 
bands (left panel). North is up and east to the left. The corresponding z'-image (right panel) illustrates the 2-D distribution of the 44 
sources in the GROND photometric sample with a photo-z that is consistent with the spectroscopic redshift, i.e., at 1 .01 < z < 1 .23 
(empty circles). These 44 photometric members are colour-coded according to their spectro-photometric classification as E (red), 
Sbc (yellow), Scd (green), Im (cyan), SB (blue) or QSO (purple). The three spectroscopic members among them are marked with 
small, black empty circles. Two concentric black circles with radii of 6" and 60" mark the X-ray position and bona-fide region 
of XMMU J0338. 7+0030, respectively. Magenta contours represent the weak, marginally extended X-ray emission of the source 
determined in Sect. 2.1.2. This emission is consistent with the projected distribution of the photometric cluster members. However, 
the candidate BCG (thick red empty circle) is offset by ~ 43" towards the east-south-east. 

total magnitude in the z' band but computed total magnitudes 
in the g',r',i', J, H, Ks bands from the corresponding 3" aper- 
ture magnitudes, corrected to the Kron aperture. This second set 
of total magnitudes yielded photometric redshifts that are better 
consistent with the spectroscopic ones for the three confirmed 
members of XMMU J0338.7+0030 with GROND photometry 
(Sect. 3.2.1). At the same time the distribution of these three 
spectroscopic members in the i' - z' vs z' colour-magnitude di- 
agram is much better consistent with the red-sequence of the 
cluster RDCSJ09 10+5422 at z = 1.106 (Mei et al. 2006a), as 
discussed in Sect. 3.2.2. For these reasons, we will make use of 
the aperture-corrected magnitudes in the g', r', i', J, H, Ks bands 
and the Kron magnitude in the z' band for all sources (method 
B), and refer to them as total magnitudes hereafter. 

Figure 7 illustrates the relation between the computed to- 
tal magnitude and its uncertainty in each of the GROND 
g', r', i', z', J, H, Ks channels. The regular behaviour observed for 
the GROND g', r', i', J, H, Ks broad bands is due to the fact that 
two sources hardly exhibit exactly the same measured value for 
their aperture fluxes. For the GROND g', r', i', J, H, Ks channels 
the uncertainty on the total magnitude is dominated by the uncer- 
tanty on the aperture correction at bright magnitudes and by the 
measurement error at faint magnitudes. The transition between 
these two regimes corresponds to a total magnitude that depends 
on the filter. 



in the z'-band Kron aperture for all bands. These aperture cor- 
rections (in the standard SExtractor format) amount to 

g ' = g '(3") - 0.1 143(+0.0094) , 

r' = r'(3")-0.1366(±0.0046), 

., _ /'(3»)_o.l219(±0.0051), 

z ' = z '(3")-0.1254(±0.0040), 

/ = 7(3") - 0.1613(±0.0175) , 

H = H(3") - 0.1510(±0.0197) , 

Ks = Ks(3") - 0.1463(±0.0307). 

We have then adopted two operational definitions of to- 
tal magnitudes and investigated their performances given the 
characteristics of the available imaging. As a first prescription 
(method A), we adopted the Kron magnitude in the z' band as 
the reference total magnitude and computed total magnitudes 
in the g',r',i', J, H, Ks bands from the reference total magni- 
tude plus the corresponding colour term measured within a 3" 
aperture and corrected to the Kron aperture. As a second pre- 
scription (method B), we adopted the Kron magnitude as the 
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Fig. 5. Comparison of the Petrosian magnitudes of stars (solid circles) and galaxies (empty circles) with GROND and SDSS photom- 
etry. The photometry in the GROND g', r', i', z' bands is calibrated with that of the stars in common with SDSS in the corresponding 
SDSS bands (panels a, b, c, d, respectively). Saturated stars and stars with poor photometry were flagged out. The asterisks indicate 
that the reference magnitudes are expressed in the GROND filter system. In each panel the red solid line represents equality. 
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Fig. 6. Comparison of the Petrosian magnitudes of stars (solid circles) after their photometry in the GROND J, H, Ks bands is 
calibrated with that in the corresponding 2MASS bands (panels a, b, c, respectively). Saturated stars and stars with poor photometry 
were flagged out. The asterisks indicate that the reference magnitudes are expressed in the GROND filter system. In each panel the 
red solid line represents equality. 



2.3.3. The photometric catalogue: star/galaxy separation 

As an alternative to the previous selection criteria, galaxies (ex- 
tended sources) can be separated from bona-fide stars (point-like 
sources) on the basis of their g' -z' and z' -Ks colours, in analogy 
with the BzK criterion (Daddi et al. 120041 cf. e.g. Galametz et al. 
2009). We adopted the formulation of the colour criterion given 
by Galametz et al. (their fig. 5) and applied the empirical colour 
transformation between the SDSS photometry and the Johnson- 
Cousins photometric system of Jordi, Grebel & Ammon (2006) 

g-B = 0.370 x (B - V) - 0.124 

plus suitable corrections for the different filter transmission func- 
tions and magnitude systems. The ensuing threshold correspond- 



ing to the GROND g'z'Ks criterion (short-dashed line in Fig. 8) 
was contrasted against the behaviour of our catalogue sources 
classified as stars on the basis of the SExtractor stellarity in- 
dex in Sect. 2.3.2. This comparison prompted the introduction 
of a shift equal to +0.27 mag in z' - Ks (solid line in Fig. 8). 
Therefore, the g'z'Ks criterion hereafter adopted to single out 
bona-fide stars is 

z - Ks (3") < 0.3 x [g' - z' (3")] - 0.1, (1) 

where we adopted 3"-aperture measurements. Comparison with 
the source classification provided by the use of the photometric 
redshift technique (Fig. 13a in Sect. 3.2.1) reveals the presence 
of a few extragalactic point-like sources with g' -z' (3") ~ 1 and 
z' - Ks (3") ~ that are misclassified as stars in Fig. 8. 
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Fig. 7. Uncertainty of the computed total photometry in each of the g', r', i', z', J, H, Ks broad-bands of GROND as a function of the 
total magnitude of a source (panels a, b, c, d, e, f, g, respectively). In each panel the short-dashed, solid and long-dashed red lines 
represent the lcr, 3<x and 5<x limits, respectively. Only sources brighter than the limiting magnitude in a given band are reproduced. 
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As an independent test of the overall goodness of the pre- 
vious criterion, in Fig. 9, we compared the distribution of the 
bona-fide stars identified through Eq. 1 in the / — Ks vs z' - Ks 
colour-colour diagram with the loci occupied by spectrally clas- 
sified stars with SDSS and 2MASS photometry (Galametz et al. 
2009, their fig. 1) after introducing suitable corrections for the 
different filter transmission functions. As discussed by Galametz 
et al., most of the stars in a matched 2MASS/SDSS catalogue 
populate either the region of the J - Ks vs z' - Ks colour-colour 
diagram above the short-dashed line and tend to be redder than 
zf — Ks = if their individual spectral types are later than K5, 
or the region below the short-dashed line and across/to the right 
of the solid line if they are earlier than K5. The consistent be- 
haviour we find confirms the validity of Eq. 1 . 



2.4. Optical spectroscopy with FORS2 



2 

(3") [AB i 



Fig. 8. Distribution of all sources with fluxes in the GROND 
g',z',Ks channels above the respective l<x thresholds in the 
z' - Ks vs g' - z' colour-colour diagram. Magnitudes within 
a 3"-wide circular aperture are used. Sources identified as point- 
like on the basis of a SExtractor stellarity index larger than 
0.98 are marked with red asterisks. Conversely, extended sources 
within (beyond) 1' from the position of XMMU J0338.7+0030 
are marked with filled (empty) circles. The short-dashed and 
solid lines represent, respectively, the computed and adopted 
thresholds to select bona-fide stars (see text). 
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Fig. 9. Distribution of all bona-fide stars, selected according to 
Eq. 1, in the J - Ks vs z' - Ks colour-colour diagram, where 
3"-aperture photometry is used. This distribution is consistent 
with the loci occupied by spectrally classified stars with SDSS 
and 2MASS photometry (cf. fig. 1 in Galametz et al. 120091 ) . As 
shown by Galametz et al., stars with a spectral type later than 
K5 populate the region above the short-dashed line and tend to 
exhibit z' - Ks > 0. Conversely, stars with a spectral type earlier 
than K5 populate the region blueward of the horizontal short- 
dashed line and define a linear colour-colour relation (solid line) 
with extensions towards very red z' - Ks colours. 



2.4.1. Observations and data reduction 

The SDSS DR7 archive contains no spectroscopy for the re- 
gion of the sky encompassing XMMU J0338. 7+0030. However, 
sparse spectroscopic information is available from XDCP, 
since the neighboring distant cluster XMMU J0338. 5+0029 
(Fassbender et al. 2011c) was observed on November 7, 2007 
(programme ID 079.A-0634, service mode) for 3 hr with 
FORS2 (Appenzeller et al. 1998). This is a multi-mode (imag- 
ing, polarimetry, long slit, MOS) instrument working at wave- 
lengths of 3300-11000 A, mounted on the Cassegrain focus of 
VET UT1 in Paranal, Chile. 

A single mask with 32 slits of 1 "-width each was designed 
and used in the FORS2 MXU-mode (Mask eXchange Unit) with 
the standard resolution collimator, which yields a 6.8' x 6.8' FoV 
with a scale of 0.25"/pixel. Coupled to the 3001 grism and the 
OG590 order blocking filter, this provided a resolution R = 660 
and a wavelength coverage of 6000-10000 A. Six single expo- 
sures were taken during the night for the same set-up, yielding 
a total net integration time of 8400 s; the seeing ranged between 
1.2" and 1.3" (FWHM). Because XMMU J0338.7+0030 was 
the secondary target of these FORS2 observations, only 30% of 
the slits were used to probe galaxies in its region. As primary 
and secondary spectroscopic targets were chosen, respectively, 
very red and bluer objects in the z - H vs H colour-magnitude 
diagram built out of the OMEGA2000 photometry (Fig. 3). As a 
result, only five slits out of a total of 32 turned out to correspond 
to objects in the region of XMMU J0338.7+0030 later imaged 
with GROND, for which photometry was extracted. 

A new version of the VIMOS Interactive Pipeline and 
Graphical Interface (VIPGI; Scodeggio et al. 2005), modified 
to make it usable on FORS2 multi-slit spectra, was used. A de- 
tailed discussion of this new multi-spectra reduction pipeline is 
deferred to a forthcoming paper (Nastasi et al. in preparation). In 
brief, the stacked 2-D raw data were flat-fielded, bias-subtracted 
and calibrated in wavelength by means of a helium-argon lamp 
used as a reference line spectrum. The final dispersion solu- 
tions had a root-mean-squares calibration uncertainty lower than 
1 A/pixel. The continuum of each spectrum was corrected for 
the sensitivity function of the FORS2 instrument, and the 1-D 
reduced spectrum was extracted from each slit for the redshift 
determination. Each 1-D reduced spectrum was both visually in- 
spected by means of the graphical tools available in VIPGI and 
cross-correlated with template galaxy spectra using the IRAF 
package RVSAO (Kurtz & Mink 1998). 
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Fig. 10. Histogram of the measured galaxy spectroscopic red- 
shifts in bins of Az = 0.01. An overdensity at z — 1.1 is evident. 
Only five of these redshifts belong to objects in the GROND 
photometric catalogue. 



2.4.2. The spectroscopic redshift of XMMU J0338.7+0030 

Figure 10 shows the redshift distribution of the galaxies at 
z > 0.5 that were observed with FORS2 in a 6.8' x 6.8' re- 
gion of the sky that contains both XMMU J03 3 8. 5 +0029 (pri- 
mary target) and XMMU J0338.7+0030 and that are not mem- 
bers of XMMUJ0338.5+0029 (see Fassbender et al. 2011c). 
Four galaxies exhibit consistent redshifts and are therefore iden- 
tified as spectroscopic cluster members. As a consequence, 
XMMU J0338. 7+0030 is at a mean spectroscopic redshift < 
z >= 1.097 ± 0.002 (lcr, where the uncertainty is determined 
from a boostrap analysis combined with the biweight location 
estimator. These four spectroscopic members lie within 2.3' 
from the X-ray position of their host cluster, which corresponds 
to a projected distance of ~ 1.1 Mpc for the assumed cosmol- 
ogy (Sect. 1). Their spectra are shown in Fig. 11, where the most 
prominent spectral features are labelled; their main properties 
are listed in Table 1. 

Out of these four spectroscopic members, three (with ID=2, 
ID=15 and ID=16) have counterparts in the GROND photomet- 
ric catalogue and lie within 1' from the original X-ray position 
of XMMU J0338.7+0030 (see Fig. 4). The object with ID=2 can 
be identified as a passively evolving galaxy (likely an early-type 
galaxy) on the basis of its spectral properties. Conversely, the 
spectroscopic member with ID=15 exhibits robust evidence of 
on-going star-formation activity, as traced by the (unresolved) 
[O II] doublet at A = 3727 A. This activity is not as vigorous as 
in the spectroscopic member with ID=12 (modulo dust attenua- 
tion) though. The object with ID=16 has the least robust deter- 
mination of photo-z and a controversial evidence of [O II] line 
emission. It is located in correspondence of the X-ray centroid 
of XMMU J0338.7+0030 determined in Sect. 2.1.1. 

The absence/presence of the [Oil] line and its equivalent 
width (EW) are consistent with the measured values of the 



Fig. 11. Spectra of the four cluster members smoothed with a 
seven-pixel boxcar filter. Black, dashed lines indicate the ex- 
pected positions of prominent spectral features at z = 1.1. 
Features caused by telluric absorption (top panel) or sky emis- 
sion (bottom panel) are plotted in red. 

strength of the 4000 A break (D4000) as shown in Table 1. In 
particular, a value of D4000 ~ 1 .9 is consistent with a minimum 
light-weighted age of the stellar populations of ~ 3 Gyr, accord- 
ing to the set of models with different star-formation histories 
and metallicities in Gallazzi et al. (2005) (A. Gallazzi private 
communication). Together with the lack of detected [Oil] line 
emission, this confirms the spectroscopic member with ID=2 
as a galaxy dominated by old, passively evolving stellar pop- 
ulations (i.e., an early-type galaxy). Minimum light-weighted 
ages of the stellar populations equal to ~ 1 Gyr and ~ 1 .4 Gyr 
seem to characterize the spectroscopic members with ID=15 and 
ID =16, respectively. These galaxy are likely evolved systems 
with ongoing star-formation activity. Finally, the spectroscopic 
member with ID=12 exhibits not only the strongest ongoing 
star-formation activity but also, consistently, the lowest value of 
D4000, which sets the minimum light-weighted age of its stellar 
populations to be equal to ~ 0.3 Gyr. 

3. Results 

3. 1. Galaxy number counts 

In order to illustrate depth and completeness of the photometric 
catalogue extracted from the GROND optical/near-IR imaging 
(Sect. 2.3.3), the number of galaxies per 0.5 mag-bin and unit 
area for each of the i',z',Ks bands is reproduced in Fig. 120 
Galaxies are selected as in Sect. 2.3.3. Statistical uncertainties 
are determined from the approximate formulae for confidence 
limits based on the Poisson and binomial statistics presented in 
Gehrels ( [19861 

6 Corrections for masked areas are introduced here; they correspond 
to 0.8% of the total area, which is equal to 4.64 x 10~ 3 deg 2 . Since these 
corrections are negligible, they will not be introduced later on. 
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Table 1. Summary of the main properties of the four galaxies observed with OMEGA2000 and FORS2 identified as spectroscopic 
cluster members. 



ID 


RA (J2000) 


Dec (J2000) 


H 


z-H 




z 


D4000 


[Oil] EW 




hh mm ss 


dd mm ss 


Vega mag 


Vega mag 








A 


2 


03 38 44.136 


+00 30 29.95 


19.39 


2.84 


1.09437 


± 0.00090 


1.901 +0.052 




12 


03 38 43.026 


+00 32 18.55 


19.03 


1.98 


1.09787 


± 0.00040 


1.194 + 0.036 


-32.908 ± 1.380 


15 


03 38 40.912 


+00 30 10.66 


18.96 


2.27 


1.09641 


± 0.00027 


1.633 + 0.028 


-7.952 ± 0.962 


16 


03 38 44.299 


+00 30 03.71 


20.29 


3.20 


1.09930 


± 0.00044 


1.690 + 0.131 


-9.291 ±3.873 



The GROND i'-, z'- and Ks-band galaxy number counts 
are compared with those determined from observations of deep 
fields or survey areas like the Hawaii Hubble Deep Field North 
(Hawaii HDF-N, Suprime-Cam I- and z'-band - Capak et al. 
120041 ). the Cosmic Evolution Survey (COSMOS, ACS WFC 
F814W filte r - Le authaud et al. 120071 Suprime-Cam i + band - 
Capak et al. l2007l WIRCam Ks band - McCracken et al. EOlOl 
the Calar Alto Deep Imaging Survey (CADIS, OMEGA Prime 
K' band - Huang et al. |20"0TT ), t he Sub aru Deep Field (SDF, 
CISCO K' band - Maihara et al. 1200101) , the infrared comple- 
ment of the VI MOS-V LT deep survey (VIRMOS, SOFI Ks band 
- Iovino et al. 120531) . the FLAMINGOS Extragalactic Survey 
(FLAMEX, FLAMINGOS Ks band - Elston et al. l2006l and the 
UKIRT Infrared Deep Sky Survey Ul tra-dee p survey (UKIDSS 
UDS, UKIRT Ks band - Hartley et al. l2008b . 

The present GROND z'-selected sample of galaxies has a 
50% completeness limit equal to z' compl ~ 25.3, as determined 
from the extrapolation of the galaxy number counts in Capak et 
al. (|2004). This completeness magnitude corresponds to a ~ 2<x 
fiux detection (cf. Fig. 7d). The 50% completeness limit is equal 
to ~ 24.8 AB mag in i'-band, which corresponds to a ~ 5<x flux 
detection. Finally, the GROND photometric sample is 50% com- 
plete down to Ks compl ~ 22.9, which corresponds to a ~ l<x 
flux detection. This completeness magnitude also corresponds 
to ~ Ks* + 2.4, where Ks* = 20. 5^ 4 is the characteristic mag- 
nitude (Schechter 1976 ) of the Ks-band luminosity function of 
cluster galaxies at z ~ 1.2 (Strazzullo et al. 2006). 

Comparison of the different galaxy number counts in Fig. 12 
shows an overall agreement between our results and those in the 
literature down to the completeness limit of the GROND photo- 
metric sample of galaxies in each band. This is important when 
considering the differences in the sky area that is probed by dif- 
ferent observations (leading to cosmic variance), the operational 
definition of total magnitude and the adopted galaxy /star sepa- 
ration criteria. These differences mostly affect the bright end of 
the galaxy number counts, which corresponds to magnitudes ap- 
proximately brighter than 20 AB mag in Fig. 12. Furthermore, 
our galaxy number counts refer to the same sample extracted in 
the z'-band, whereas the galaxy number counts in the literature 
refer to individual samples extracted in the i'-, z'- or K-band. 

In spite of the overall agreement, our number counts exhibit 
an excess (by a factor of almost two) in the magnitude range be- 
tween Ks comp] - 2 and ^""P 1 - 0.5 (i.e., 20.8 < Ks < 22.3) 
with respect to all others. This magnitude range corresponds to 
~ Ks* + 0.3 < Ks < K s* + 1.8 for cluster galaxies at z ~ 1.2 
(Strazzullo et al. 2006). Consistently, a slight excess across a 
broader magnitude range (i.e., approximately at magnitudes be- 
tween 21.5 and 23.5 AB mag) is visible from comparison of the 
galaxy number counts for the z' band. No clear excess is appar- 
ent in the galaxy number counts for the i'-band. This behaviour 
is interpreted as an indication for the presence of an over-density 

7 We plot the completeness-corrected, K'-band galaxy number counts 
provided by Maihara et al. ( 12001 \ . 



of galaxies associated with the weak, marginally extended X-ray 
emission corresponding to XMMU J0338.7+0030. 

3.2. Photometric-redshifts 

The photometric redshift technique enables us to assess the 
membership to XMMU J0338.7+0030 for the 832 sources in the 
GROND photometric sample. Distances to individual sources 
were estimated through the publicly available package for 
PHotometric Analysis for Redshift Estimations, le Phar^ (S. 
Arnouts & O. Ilbert), which is a set of Fortran programmes to 
compute photo-z's using the standard SED fitting technique. 

We selected the large pool of observed SEDs of stars (254) 
and the mixture of observed and synthetic SEDs for QSOs (28) 
that are provided by le Phare, as it is standard. However, we ini- 
tially adopted two sets of templates to describe observed SEDs 
of galaxies and investigated their performances. 

The first set was used for the COSMOS photo-z paper (Ilbert 
et al. 120091 ). It contains 31 templates: the SEDs numbered from 
1 to 7 correspond to galaxies that can be morphologically classi- 
fied as elliptical (E), the SED number 8 corresponds to a lentic- 
ular (or SO) galaxy, SEDs numbered from 9 to 12 correspond to 
bulge-dominated disc galaxies (i.e, Sa and Sb), SEDs numbered 
from 13 to 19 correspond to bulge-less disc galaxies (i.e, Sc, Sd, 
and Sdm), SEDs numbered from 20 to 3 1 describe starburst (SB) 
galaxies. We followed the indicated prescriptions for the model- 
dependent attenuation^ of galaxy templates, so that a "Prevot ex- 
tinction law" (Prevot et al. 1984) was applied to SEDs numbered 
from 13 to 23, a so-called "Calzetti attenuation law" (Calzetti et 
al. 2000) and two "modified" Calzetti attenuation laws (see Ilbert 
et al. 120091 ) were applied to SEDs numbered from 23 to 31. This 
choice accounts for the fact that high-z star-forming galaxies are 
characterized by extinction laws that differ at rest-frame ultravi- 
olet wavelengths, as it is in the Local Group (Noll et al. 2009 
and references therein). In all cases, the reddening E(B - V) was 
enabled to range between and 0.5 mag in steps of 0.05 mag. 

The second set of galaxy templates corresponds to the 66 
SEDs that are used in Ilbert et al. (2006) for the analysis of the 
CFHTLS Deep fields. In the CFHTLS set, the ratio between the 
numbers of E and SB templates is reversed with respect to the 
COSMOS set: E, Sbc, Scd, Im and SB galaxies are described 
by 22, 17, 12, 1 1, and 4 SEDs, respectively. For Scd, Im and SB 
galaxy templates, a Prevot et al. (1984) extinction law is adopted, 
with values of E(B-V) ranging from to 0.25 with a step of 0.05 
consistently with the indications in the le Phare documentation. 



8 www.oamp.fr/people/arnouts/LE_PHARE.html 

9 We note that the extinction curve describes the combined absorp- 
tion and out-of-the beam scattering properties of a mixture of dust 
grains of given size distribution and chemical composition in a screen 
geometry as a function of wavelength; the attenuation function is the 
combination of the extinction curve with the geometry of a dusty stellar 
system, in which a substantial fraction of the scattered light is returned 
to the line of sight (e.g. Pierini et al. 2004b). 
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Fig. 12. Comparison of the galaxy number counts determined 
in the region of XMMU J0338.7+0030 imaged with GROND 
for the i' (top), z' (centre) and Ks band (bottom) and the cor- 
responding ones from deep fields or survey areas. Error bars are 
shown only for our points. In each panel, the short-dashed, solid 
and long-dashed lines represent the 1 <x, 3 <x and 5 <x thresh- 
olds of our photometry, respectively. The GROND galaxy num- 
ber counts exhibit a clear excess for 20.8 < Ks < 22.3. 
14 



For each library of galaxy SEDs, we considered two sepa- 
rate families of models, by including the contribution to the total 
flux from line emission or not. We opted to do this because there 
is independent evidence that the bright end of the galaxy lumi- 
nosity function of a high-z cluster can be populated by [Oil] 
line emitters (as for XMMU J0302.2-0001 atz = 1.186, Suhada 
et al. 2011), sometimes with vigorous starburst activity (as for 
XMMU J1007.4+1237 at z = 1.555, Fassbender et al. 201 lb). 

Finally, for galaxy templates as well as for QSO templates, 
we considered the redshift range 0-6 with a step of 0.02 in z and 
the range -24 < M z > < -8 in the z'-band absolute magnitude. 

For all stellar, QSO and galaxy templates, synthetic photom- 
etry was computed for all GROND bands in the AB magnitude 
system and confronted with the observed photometry. Both sys- 
tematic (e.g., caused by calibration uncertainties) and statisti- 
cal (i.e., caused by measurement uncertainties) errors were taken 
into account, the latter being magnified by 50% (in mag). 

The photometric redshift code le Phare, based on a simple^ 2 
fitting method, yields the best-fitting template in each class, the 
best-fit value of z, the 68% range of the photo-z solutions and the 
probability per redshift bin, in particular. All ouputs for individ- 
ual detected sources were inspected by eye. A source was classi- 
fied as a star if the best-fitting stellar template gave the lowest^ 2 
value. Conversely, it was classified as a QSO if the best-fitting 
QSO template gave a^ 2 value that was better by a factor of two 
than those given by the best-fitting stellar and galaxy templates. 

As a first step, we checked the reliability of the publicly 
available transmission functions of the GROND dichroic filters 
uploaded into le Phare by comparing the distributions of the 
library stars and the stars observed by us with the most accu- 
rate photometry in a large suite of colour-colour diagrams (see 
Fig. 13 for examples). Eventual problems affecting the GROND 
photometry of each star in individual bands (e.g., owing to poor 
S/N and deblending) were inspected: if there were compro- 
mised photometric data, they were excluded. As a result, we 
found an overall good agreement between observed and syn- 
thetic stellar tracks, which required the addition of very small 
shifts (-0.05 mag, +0.025 mag, +0.025 mag) to the synthetic 
photometry in the r',i',H bands, respectively. Furthermore, the 
colour criterion adopted to statistically select stars from g', z', Ks 
photometry only (Eq. 1) was nicely confirmed (cf. Figs. 8 and 
13a). 

As a second step, we computed photometric redshifts for 
two photometric catalogues with the same sources as part of the 
test to determine the best operational definition of total mag- 
nitude for our data. This test was executed with the COSMOS 
set of galaxy templates excluding the line emission contribu- 
tion; in particular, it addressed the five galaxies with GROND 
photometry and a spectroscopic redshift (Sect. 2.4.1). Overall 
the uncertainty on the ensuing two sets of photometric redhifts 
obtained for these five galaxies is similar, as well as the best- 
fit solutions yielded for the two foreground galaxies. However, 
a mean photometric redshift equal to ~ 1.22 + 0.25 (lcr) or 
~ 1.02 ± 0.15(lcr) is obtained for XMMU J0338.7+0030 when 
total magnitudes were computed using the operational defini- 
tions A and B in Sect. 2.3.2, respectively. Because this cluster 
is spectroscopically confirmed to be at z = 1.1 (Sect. 2.4.2), the 
operational definition of total magnitude B has to be favoured. 
This choice is also supported by the distribution of the three 
spectroscopic members in the i' - z' vs z' colour-magnitude dia- 
gram (Fig. 15). Therefore, we only discuss results from the fidu- 
cial photometric catalogue where total magnitudes are computed 
with method B hereafter. 
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Fig. 13. Five colour-colour diagrams for stars in the le Phare library (black filled circles) and the 15 to 23 objects detected at all 
GROND bands, with the most accurate photometry in the relevant bands and robustly classified as stars by le Phare using the full 
seven-channel photometry (red asterisks). These observed stars nicely fall onto the synthetic stellar sequences, after the synthetic 
photometry is slighlty corrected by -0.05, +0.025 and +0.025 mag in the r', i' and H bands, respectively. We recall that photometric 
redshifts are particularly sensitive to the 4000 A break in the spectrum of a galaxy, which spans the entire wavelength range of thg 
z' broad-band filter for z = 1 . 1-1 .6. 
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Finally, we compared the performances of the COSMOS and 
CFHTLS sets, whether the emission line contribution to the tem- 
plates of star-forming galaxies is included or not. We based our 
judgement on the comparison of the best-fit solutions and their 
uncertainties but also on the capability of the set of choice to 
recover the redshifts of the five galaxies with available spec- 
troscopy (see Sect. 2.4.1). This figure corresponds to only ~ 
0.6% of the GROND photometric sample of 832 sources; the en- 
suing span in redshift is very narrow: z — 0.85-1 . 1 . Nevertheless, 
the results obtained for the three spectroscopic members of 
XMMUJ0338.7+0030 (see Sect. 2.4.2) gave us useful indica- 
tions. The use of the COSMOS galaxy templates yielded a mean 
photometric redshift equal to ~ 1.02 + 0.15 for this cluster, 
whether emission lines were included or not. Conversely, the use 
of the CFHTLS galaxy templates yielded values of ~ 1.14 ±0.1 1 
and ~ 1.12 + 0.09, respectively. Given the better agreement with 
the spectroscopic redshift of XMMU J0338.7+0030 (< z >= 
1 .097 + 0.002) and the lower uncertainty, hereafter we discuss 
results obtained from the use of the CFHTLS set of templates 
without emission line contributior{3 

As a consequence of this choice, out of the 832 entries in the 
GROND photometric catalogue, 174 sources do not have a (re- 
liable) photo-z estimate (owing to very poor photometry or mul- 
tiple peaks in the photo-z probability distribution function), 64 
sources are identified as stars, 158 sources exhibit SEDs that are 
best-fit by a QSO template, 436 sources are identified as galax- 
ies. On one hand, the past large percentage of objects identified 
as QSOs reflects our arbitrary classification criterion of the stel- 
lar, galaxy and QSO templates that best-fit the SED of a given 
object. On the other hand, it reflects the reduced capability of 
discriminating between best-fit templates for the numerous ob- 
jects with incomplete photometric information. For the majority 
of the objects identified as QSOs, the photometric redshift is re- 
assuringly consistent, if not similar, to the value associated with 
the best-fit galaxy template. 

The quality of the best-fit photometric redshifts obtained 
with CFHTLS galaxy templates not including emission lines is 
discussed in Appendix A. The distribution of these photo-z solu- 
tions is often single peaked and narrow when accurate photome- 
try is available for all seven GROND channels. 

3.3. Optical identification of the high-z cluster 

3.3.1. The photometric-redshift distribution function 

We computed the surface density of extragalactic sources within 
or beyond a cluster-centric distance of 1', as well as across the 
entire 3.9 x 4.3 arcmin 2 region of the sky imaged with GROND, 
as a function of the photometric redshift (Fig. 14a). A simple vi- 
sual inspection reveals that for the redshift bins 1.0-1.1, 1.2-1.3 
and 1.5-1.6, the density of sources within the bona fide cluster 
region exceeds that computed in the adjacent field by a factor of 
2. An equivalent over-density of V - /Ti-selected extremely red 
objects (EROs£3 exists in each of the last two bins (Fig. 14b). 



10 This suggests that the COSMOS galaxy templates are less ad- 
equate to represent the characteristics of the galaxy populations of 
XMMU J0338. 7+0030. An analogous point was made by Guennou et 
al. (2010), who discussed photometric redshifts along lines of sight to 
10 clusters at spectroscopic redshifts between 0.48 and 0.79 obtained 
with the same COSMOS galaxy templates and photo-z code used here. 

11 Colours redder than I - K = 4 Vega mag statistically select EROs 
at z ~ 1-2 (Pozzetti & Mannucci 2000 and references therein; see also 
Pierini et al. 2004a 2005 ). For the GROND channels, this corresponds 
to i' - Ks — 2.8 in the AB magnitude system. 



candidate cluster region 
field 

entire FoV 
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field EROs 
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Fig. 14. Number density as a function of the photometric redshift 
for the extragalactic sources with a photo-z less than 2 in the en- 
tire 3.9 x 4.3 arcmin 2 region of the sky imaged with GROND 
(black solid line), for those among them inside the bona fide 
cluster area (red solid line) or in the surrounding field (blue 
short-dashed line) - panel a. The analogous distribution func- 
tions for the sub-sample of EROs, expected to be at z ~ 1-2 (see 
text), are also reproduced (panel b). For the redshift bins 1.0- 
1.1, 1.2-1.3 and 1.5-1.6, the density of sources within a cluster- 
centric distance of 1' exceeds that computed in the adjacent field 
by a factor of 2. In correspondence of the last two redshift bins, 
there is a clear excess of EROs. 



These EROs are found in the redshift range where they are ex- 
pected to be, which gives additional confidence in the quality of 
our photo-z's, at least in a statistical sense. 

In support of this conclusion, the spectroscopic members 
with ID=2, ID=15 and ID=16 have photometric redhift{3 
equal to 1.01 10' ^, 1.1711 j j^ and 1.1820j ^, respectively. 
Furthermore, the morphological type of the best-fitting template 



12 Here we give the best-fit photo-z together with the values that 
bracket the 68% range of the photo-z solutions. 
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Fig. 15. i'—zf vs z' colour-magnitude diagram for all GROND sources in the XMMU J0338.7+0030region with flux detections in the 
i', z' bands at a significance level > lcr, whatever the classification of the source by Le Phare (i.e., stars are included). Photometric 
errors (lcr) are also shown. Candidate members of XMMU J0338. 7+0030 are reproduced with red (green) filled circles if they 
are within (outside) the bona fide cluster region (panel a). The three spectroscopic cluster members in the GROND photometric 
catalogue are marked with blue filled circles there. Green asterisks represent the two QSO at a cluster-centric distance greater 
than 1'. Alternatively (panel b), photometric cluster members classified as galaxies are colour-coded according to their spectro- 
photometric type, which is E (red), Sbc(yellow), Scd (green), Im (cyan) or SB (blue). In each panel, the red short-dashed, solid and 
long-dashed lines represent the lcr, 3cr, and 5cr flux thresholds of our photometry (cf. Sect. 2.3.2). The red dotted-short-dashed line 
represents the completeness limits in the i',z' bands, as obtained from Fig. 12. As a reference, the blue solid line reproduces the 
red-sequence defined by the inner regions of the 31 E+SO+SO/a galaxies that are located within 1' from the X-ray centroid of the 
cluster RDCS J09 10+5422 at z = 1.106 (Mei et al. 2006a). 



corresponds to an E galaxy for the spectroscopic members with 
ID=2 and ID=16 and to an Scd galaxy for the spectroscopic 
member with ID=15. This spectro-photometrically determined 
classification is nicely consistent with the spectroscopic classi- 
fication in Sect. 2.4.2. Furthermore, it is reasonable to conclude 
that the two peaks atz~ 1.0-1.1 andz ~ 1.2-1.3 in Fig. 14a are 
caused by true cluster members, at least in part. 

Therefore we selected as photometric (also candidate or fidu- 
cial hereafter) cluster members those sources with a photo-z in 
the range 1.0110-1.2317, which is centred on the mean pho- 
tometric redshift of XMMU J0338.7+0030 (< z >= 1.1214 ± 
0.0957, lcr), contains the previous three spectroscopic members 



and, thus, has a width equal to ±1.15cr (or ±0.1 < z >). Forty- 
four sources satisfy our criterion: fourteen (36) are within 1' (2') 
from the X-ray position of XMMU J0338.7+0030. Only two of 
these 44 candidate cluster members are identified as QSOs ac- 
cording to their best-fit templates; they are at cluster-centric dis- 
tances greater than 1'. The remaining 42 fiducial cluster mem- 
bers are spectro-photometrically classified as E (7), Sbc (6), Scd 
(11), Im (8) or SB (10) according to the templates that best fit 
their individual SEDs and provide their fiducial photometric red- 
shifts. As discussed in Appendix A, the photometric redshifts as- 
sociated with these 42 galaxies are, in general, well-behaved and 
exhibit low uncertainties. 
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3.3.2. The i' - z' vs z' colour-magnitude diagram: evidence 
for a red sequence 

We now investigate if the previously selected candidate clus- 
ter members consistently define a red sequence in a suitable 
colour-magnitude diagram. We considered the i',z' bands be- 
cause the rest-frame 4000 A break spans the entire wavelength 
range of the z' band as the redshift of a source increases from 
1.1 to 1.6. Figure 15 shows the distribution of all sources of 
the GROND photometric sample in the i' - z' vs. z' colour- 
magnitude diagram. In the upper panel, the photometric mem- 
bers of XMMU J0338.7+0030 are marked in red or green if they 
fall inside or outside the bona-fide cluster region, respectively. 
In the lower panel, the same sources are colour-coded accord- 
ing to the best-fit template (if they are identified as galaxies). 
About a dozen candidate cluster members identified as galax- 
ies (including the three spectroscopic cluster members) sug- 
gest the existence of a red locus in the i' - z' vs. z' colour- 
magnitude diagram: i' -z' slightly decreases from ~ 1.1 AB mag 
to ~ 1.0 AB mag as z' increases from ~ 21.5 AB mag to 
~ 24.0 AB mag. Furthermore, half of these galaxies are iden- 
tified as ellipticals. This is consistent with the existence of a red 
sequence (Arimoto & Yoshii 1987). 

As a reference, the red sequence of RDCS J09 10+5422 at 
z = 1.106 (Mei et al. 2006a) is reproduced in both panels of 
Fig. 15. Mei et al. obtained high-resolution imaging of this mas- 
sive cluster with the Hubble Space Telescope Advanced Camera 
for Surveys (HST ACS) in the F775W and F850LP bandpasses. 
They computed total Z850 magnitudes (SExtractor MAG_AUTO) 
and i-ns - Z850 colours within the effective radii of the individ- 
ual galaxies in order to select candidate cluster members. The 
use of these aperture colours was meant to avoid a selection bias 
introduced by the potential presence of radial colour gradients. 
Furthermore, Mei et al. fitted different (775 - zsso vs Z&50 colour- 
magnitude relations to different sub-samples of morphologically 
classified early-type galaxies. We have transformed their best-fit 
relation for the 31 E+SO+SO/a galaxies within 1' from the X-ray 
centroid of RDCS J09 10+5422 (see table 1 in Mei et al. 2006a) 
into a red sequence in the i' - z' vs z' colour-magnitude dia- 
gram by introducing corrections for the different transmission 
functions of the corresponding HST ACS and GROND filters. In 
spite of the uncertainties on these colour transformations (likely 
of the order of 0. 1 mag), the red-sequence of RDCS J09 1 0+5422 
describes the red locus of XMMU J0338.7+0030 rather well. 

However, galaxies with very different spectro-photometric 
types populate the red locus of XMMU J0338.7+0030 within 
a cluster-centric distance of 1' and not just early-type galax- 
ies as in RDCS J0910+5422. This was also observed in 
e.g. XMMUJ1230.3+1339 at z = 0.975 (see fig. 7 in 
Lerchster et al. 2011). Nevertheless the two coeval clusters 
XMMUJ0338.7+0030 and RDCS J09 10+5422 seem to host 
also analogus populations of star-forming galaxies that exhibits 
bluer i' - z' colours with respect to the red, passively evolving 
members by about 0.5 AB mag (cf. fig. 7 in Mei et al. 2006a). 
For XMMU J0338.7+0030, these galaxies constitute about half 
of the candidate members and are identified as Scd, Im or SB. 
Furthermore, those within the bona-fide cluster region span the 
magnitude range 22.9 < z' < 24.9, whereas those farther out 
can be as bright as z' ~ 21.7. These bright, star-forming, can- 
didate cluster-member galaxies might be analogous to the spec- 
troscopic member with ID=12 at a cluster-distance of 2.3' (cf. 
Table 1). Evidence of a population of bright star-forming disc 
galaxies was found in RDCS J09 10+5422 as well and within a 
cluster-centric distance of 2' (Mei et al. 2006a). 
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Fig. 16. Cumulative number density of the 44 photometric mem- 
bers of XMMUJ0338.7+0030 as a function of their cluster- 
centric distances. The mean number density of those members 
within a concentric, circular annulus delimited by radii of 1' and 
2' is reproduced together with upper and lower limits (red solid 
and short-dashed lines, respectively) computed as in Gehrels 
(1986). 



In conclusion, the photometric redshifts determined in 
Sect. 3.2 confirm that XMMU J0338.7+0030 is a cluster at z = 
1 . 1 and suggest that its galaxy population is qualitatively similar 
to that of the coeval cluster RDCS J09 1 0+5422. 



3.3.3. Projected galaxy distribution 

In the entire region of XMMU J0338.7+0030 imaged with 
GROND at all seven bands, there are 44 candidate cluster mem- 
bers, including three spectroscopic members (Sect. 3.2.1). Their 
projected distribution as a function of cluster-centric distance 
(up to 2') is reproduced in Fig. 16. There are 14 photometric 
members in the bona-fide cluster region and 22 in the concen- 
tric annulus of inner and outer radii equal to 1' and 2' that has 
a three times larger area. If one assumes this annulus to repre- 
sent the surrounding low-density region (i.e., the "field"), the 
over-density of galaxies within 1' from the X-ray position of 
XMMU J0338.7+0030 has a significance of 4.3cr. The signifi- 
cance rises to 5.3<x if one assumes that the region external to the 
circular area of 2' radius centred on the X-tay position of the 
cluster represents the coeval field environment. Considering that 
a value of R 200 equal to 88.8" (equivalent to ~ 725 kpc) can be 
estimated from the scaling relations in Fassbender et al. (201 lc), 
we conclude that the existence of an over-density of galaxies 
at the same photometric redshift of XMMU J0338.7+0030 is 
highly significant. This is consistent with the excess by a factor 
of two in the Ks-band galaxy number counts (for the entire re- 
gion imaged with GROND) at ~ 20.5-22.5 AB mag with respect 
to several computations for the field from survey areas/deep 
fields (Fig. 12 in Sect. 3.1). 

In addition, the projected distribution of the 44 photomet- 
ric members is not concentrated towards the X-ray position of 
the cluster. This is consistent only in part with the original in- 
ference from the OMEGA2000 z, H imaging that the 2-D dis- 
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tribution of the galaxies with very red z - H colours is offset 
by about 20" (see Fig. 1 and Sect. 2.2.1). Indeed, the 2-D dis- 
tribution of the 44 candidate cluster members in the z' -image 
of XMMU J0338.7+0030 (right panel of Fig. 4 in Sect. 2.3.1) 
is broadly consistent with the weak, marginally extended X-ray 
emission of the cluster. A large part of the candidate members is 
distributed in the north-south direction. This is consistent with 
the location of the spectroscopic member with ID=12, that lays 
2.3' away from the X-ray position of the cluster and to the north 
(Sect. 2.4.2). The elongation of the X-ray emission from the 
ICM suggests that XMMU J0338.7+0030 is still in an assembly 
phase. 

We identify as the brightest central galaxy (BCG) the bright- 
est photometric member, which sits on the red locus of the 
V -z' vs. z' colour-magnitude diagram of XMMU J0338.7+0030 
(Fig. 15). The BCG is offset by ~ 43", which corresponds to a 
projected distance of ~ 350 kpc from the X-ray position of the 
cluster for z — 1.1 and the adopted cosmology (Sect. 1). This 
galaxy is consistently classified as elliptical by le Phare; how- 
ever, it is very close to a bright star, which can affect its photom- 
etry. For this reason, the identified BCG was not selected as a 
target for spectroscopy with FORS2 (see Sect. 2.4.1). 

Conversely, the spectroscopic member with ID=16 is only 
2.1" away from the X-ray position of XMMU J0338.7+0030: 
this separation is smaller than the on-axis HEW (~ 14") of 
XMM-Newton. This galaxy exhibits a weak evidence for [Oil] 
line emission (cf. Fig. 1 1 and Table 1 in Sect. 2.4.2) but is iden- 
tified as an E galaxy by le Phare. With Ks = 21.19 + 0.10 
and i' - Ks = 3.36 ± 0.12 (i.e., Ks = 19.40 + 0.10 and 
i' -Ks = 4.77 + 0.1 2 in the Vega magnitude system) it also meets 
the selection criterion of EROs (footnote 11). Its classification 
as a galaxy dominated by old, passively evolving stellar popula- 
tions (i.e., an elliptical galaxy) or a dusty starburst is ambiguous 
because it exhibits J-Ks = 1.28+0.13 (i.e., J-Ks = 2.17+0.13 
in the Vega magnitude system; see Pozzetti & Mannucci 2000; 
see also fig. 13 in Pierini et al. 2004a). This spectroscopic clus- 
ter member might be associated with part of the emission de- 
tected through XMM-Newton because about 15% of the I - K- 
selected EROs with K < 20. 1 Vega mag exhibit X-ray properties 
consistent with those expected from luminous, obscured active- 
galactic-nuclei (e.g., Alexander et al. l2002l l. If so, the true X-ray 
emission centroid of XMMU J0338.7+0030 might be different. 



4. Discussion 

Sections 2.1 and 2.2 have described data reduction and analysis 
for the archival XMM-Newton X-ray observations and the near- 
IR follow-up imaging with OMEGA2000 that led to the selec- 
tion of XMMU J0338.7+0030 as a candidate cluster at z > 0.8. 
Available FORS2 observations indicate four galaxies within 2.3' 
from the X-ray position of this cluster as spectroscopic mem- 
bers. This establishes the existence of a bound system of galax- 
ies and hot plasma at < z >= 1.097 ± 0.002 (lcr) as illustrated in 
Figs. 10 and 11 and Table 1 (Sect. 2.4.2). With an estimated to- 
tal mass M 200 ~ 10 14 M , XMMU J0338.7+0030 is in between 
X-ray selected groups and clusters (Sect. 2.1.2). 

There is some tension between the spectroscopic redshift 
and the original estimate of z ~ 1.45 ± 0.15 (lcr) that was 
based on the tentative identification of red, passively evolving, 
member galaxies in the z - H vs z colour-magnitude diagram 
of XMMUJ0338.7+0030 (Fig. 3 in Sect. 2.2.2). A honest but 
wrong assumption (see Sect. 2.2.2) as well as systematics (see 
Appendix B) are at the origin of the discrepancy. 




5000 10* 1.5x10* 2x10* 2.5x10* 

A [Angstroems] 



Fig. 17. Transmission curves of the GROND g', r', i', z', J, H, Ks 
dichroic filters (black) plus normalized SEDs that represent a 
galaxy dominated by an old, passively evolving stellar popula- 
tion (red) and a very young, dusty starburst galaxy (blue). See 
the text for details. 

In addition to the spectroscopic observations with FORS2, 
we performed simultaneous optical/near-IR imaging of a re- 
gion of the sky centred on XMMU J0338.7+0030 with GROND. 
Several pieces of evidence consistently point to the existence of 
an optical counterpart to the weak, marginally extended X-ray 
emission detected with XMM-Newton in the area imaged at all 
GROND g', r', i', z', J, H, Ks bands: 

- the presence of an excess by a factor of two in the Ks- 
band galaxy number counts at ~ 20.5-22.5 AB mag with 
respect to analogous computations for the field from survey 
areas/deep fields (Fig. 12 in Sect. 3.1); 

- the consistency between the spectroscopic redshift of the 
cluster and the mean value of the photometric redshifts of the 
three spectroscopic members imaged with GROND, which is 
equal to 1.12 ±0.09 (lcr); 

- the selection of 44 sources in the photometric redshift range 
1.01-1.23 that contains the three aforementioned cluster 
members, out of which fourteen are within a cluster-centric 
distance of 1' and signpost a galaxy over-density at a signif- 
icance of 4.3<x (Sect. 3.2.1); 

- the existence of a red locus in the distribution of these 
44 photo-z selected galaxies in the i' - z' vs z' colour- 
magnitude diagram, which is consistent with the red se- 
quence of the cluster RDCS J09 10+5422 at z = 1.106 
(Fig. 15 in Sect. 3.2.2); 

- the consistency between the projected spatial distribu- 
tion of the 44 photometric cluster members and the 
weak, marginally extended emission of the X-ray source 
XMMU J0338.7+0030 (Fig. 4 in Sect. 2.3.2). 

We acknowledge that the identification of three spectro- 
scopic members out of the GROND photometric sample of 832 
sources was essential to determine the photo-z of the cluster as 
well as to select its candidate (photometric) members. Indeed, 
the photo-z distribution of the sources within 1' from the X-ray 
position of the cluster exhibits two substantial excesses at photo- 
metric redshifts of 1-1.1 and 1.2-1.3 with respect to the photo-z 



19 



D. Pierini et al.: I. The galaxy population of XMMU J0338.7+0030 at z = 1.1 




18 20 22 

Ks [AB mag] 



Fig. 18. i' - Ks vs Ks colour-magnitude diagram for all GROND sources in the XMMU J0338. 7+0030 region with flux detections 
in the i',z',Ks bands at a significance level > lcr, whatever the classification of the source by Le Phare (i.e., stars are included). 
Photometric errors (lcr) are also shown. In each panel, candidate members of XMMU J0338. 7+0030 are reproduced according 
to the corresponding panel of Fig. 15. Here the red short-dashed, solid and long-dashed lines represent the lcr, 3cr and 5cr flux 
thresholds in the i', Ks bands, respectively. The red dotted-short-dashed line represents the completeness limits in the i', Ks bands, 
as obtained from Fig. 12. The black vertical dotted and horizontal long-dashed lines reproduce the characteristic magnitude of the 
Ks-band luminosity function of cluster galaxies at z ~ 1.2 (Strazzullo et al. 2006) and the i' - Ks colour-selection threshold for 
EROs (Pozzetti & Mannucci 2000). As a reference, the blue solid and short-dashed lines reproduce the red sequences of the clusters 
RDCS J09 10+5422 at z = 1.106 (Tanaka et al. 2008) and RDCS J1252.9-2927 at z = 1.24 (Tanaka et al. |2U0"9l ), respectively. The 
candidate members of XMMU J0338. 7+0030 that are classified as elliptical galaxies (red symbols in panel b) populate a locus that 
is broadly consistent with the red sequences of the other two clusters. 



distribution of the sources populating the remaining area imaged 
by GROND at all seven bands (Fig. 14a in Sect. 3.2.1). This di- 
cotomy is caused, at least in part, by the following reasons. 

On one hand, the wavelength coverage yielded by the 
GROND dichroic filters for a z = 1.1 source is not optimal. 
This is shown in Fig. 17, where the transmission curves of the 
GROND bands are reproduced together with two synthetic SEDs 
(from Pierini et al. 2004a) redshifted to z = 1.1 and normalized 
at the effective wavelength of the z'-band. One SED corresponds 
to an SSP model with age of 4 Gyr and solar metallicity (red); it 
represents a galaxy characterized by an old, passively evolving 
stellar population. The other corresponds to a composite stellar 



population model with constant star-formation rate (SFR), age 
of 0. 1 Gyr, half-solar metallicity, attenuated by dust as for a typ- 
ical nearby starburst (blue); it represents a very young, dusty 
starburst galaxy. The redshifted 4000 A break, a fundamental 
feature in redshift determination, falls in between the i'-band 
and the z'-band. The redshifted [Oil] (A = 3727 A) emission 
line falls at the red edge of the i'-band, whereas the redshifted 
[OIII](/l = 5007 A) and Ha(A = 6563 A) emission lines fall, 
respectively, in regions of the z' and J filters with poor transmis- 
sion efficiency. 

On the other hand, some ambiguity is favoured by the photo- 
metric uncertainty yielded by the present, unoptimized GROND 
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Fig. 19. J - Ks vs Ks colour-magnitude diagram for all GROND sources in the XMMUJ0338. 7+0030 region with flux detections 
in the z', J, Ks bands at a significance level > lcr, whatever the classification of the source by Le Phare (i.e., stars are included). 
Photometric errors (lcr) are also shown. In each panel, candidate members of XMMU J0338. 7+0030 are reproduced according to the 
corresponding panel of Fig. 15. Here the red short-dashed, solid and long-dashed lines represent the lcr, 3cr and 5cr flux thresholds 
in the J, Ks bands. The black vertical dotted and horizontal dotted-long-dashed lines reproduce the characteristic magnitude of the 
Ks-band luminosity function of cluster galaxies at z ~ 1 .2 (Strazzullo et al. 2006) and the J - Ks colour-selection threshold for 
DRGs (Franx et al. 2003 see text). As a reference, the blue solid and short-dashed lines reproduce the red sequences of the clusters 
RDCS J1252.9-2927 at z = 1.24 (Lidman et al. 12004b and XMMU J2235.3-2557 at z = 1.39 (Lidman et al. 120081 respectively. 
The candidate members of XMMU J0338. 7+0030 that are classified as elliptical galaxies (red symbols in panel b) populate a locus 
that is broadly consistent with the red sequences of the other two clusters. However, candidate members that are classified as star 
forming span the same broad range of colours, especially at magnitudes fainter than ~ Ks* + 1 . 



observations, in particular at near-IR wavelengths (Fig. 7 in 
Sect. 2.3.2). It is known that the photometric redshift technique 
can be reliably extended from z ~ 1 to z ~ 2 only if deep near- 
and mid-IR photometry is available (e.g., Bolzonella, Miralles, 
Pello 2000; Ilbert et al. 2009). 

In spite of the unfortunate mapping of the rest-frame SED 
of a z - 1.1 galaxy yielded by the GROND channels and the 
limited photometric accuracy of the present near-IR data, it is 
remarkable that the photometric and spectroscopic estimates of 
the redshift of XMMU J0338.7+0030 coincide within an uncer- 
tainty of 0.09 (cf. e.g. Lerchster et al. 2011). Indeed, the avail- 
ability of spectroscopic redshifts for only 0.6% of the GROND 



sources (and in the narrow redshift range 0.85-1.1) prevented 
any meaningful training of the solutions of the photo-z code le 
Phare. 

Another exquisite possibility offered by the simultaneous 
optical/near-IR imaging with GROND is the characterization of 
the galaxy population of a high-z cluster on the basis of SEDs 
- not of a single colour-magnitude diagram - in one and the 
same study (cf. Mullis et al. 2005 and Lidman et al. 2008 for 
XMMUJ2235.3-2557; Lidman et al. 2004 and Tanaka et al. 
20 09 for RDCS J1252.9-2927; Mei et al. 2006a and Tanaka et 
al. 120081 for RDCS J0910+5422). For instance, the distribution 
in the i' - z' vs z' colour-magnitude diagram of the 44 photo- 
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z selected members of XMMU J0338. 7+0030 resembles that of 
the members of the cluster RDCS J09 10+5422 at z = 1.106, as 
determined by Mei et al. (2006a), as shown in Fig. 15. At the 
same time, the fiducial members of XMMUJ0338.7+0030 that 
are classified as elliptical galaxies populate a red locus in the 
;' - Ks vs Ks colour-magnitude diagram, which is broadly con- 
sistent with the red sequence of RDCS J09 10+5422 determined 
by Tanaka et al. (2008 ), as shown in Fig. 18b. 

The rather neat progression from red, passively evolving, el- 
liptical galaxies to blue, star-forming, late-type spirals and ir- 
regulars in Fig. 15b is not preserved in other colour-magnitude 
diagrams that are used to constrain distance and properties of 
the galaxy population of high-z clusters. The identification of a 
red-sequence and the classification of a cluster member galaxy 
as passively evolving or star forming become much more am- 
biguous in those colour-magnitude diagrams that involve broad- 
band filters that do not map the redshifted 4000 A break. This 
is demonstrated by the distribution of the fiducial members of 
XMMUJ0338.7+0030 in the i' - Ks vs Ks and J - Ks vs 
Ks colour-magnitude diagrams (Figs. 18 and 19, respectively). 
As a reference, we reproduce the red-sequences of the spectro- 
scopically co nfirme d clusters RDCS J09 10+5422 at z = 1.106 
(Tana ka eta l.|2008), RDC S J125 2.9-2927 at z = 1.24 (Lidman 
et al. 120041 Tanaka et al.l2009b and XMMU J2235.3-2557 at 
z = 1.39 (Lidman et al. l2008Fl 

A galaxy classification based on either spectro-photometry 
or morphology indicates that there are by far not as many 
early-type galaxies within a cluster-centric distance of 1' 
in XMMUJ0338.7+0030 as in RDCS J0910+5422: 3 vs 31 
(Sect. 3.2.2). A high percentage of passively evolving galaxies 
among the brightest members of a cluster should not necessar- 
ily be expected if the cluster is observed at an early stage of its 
halo assembly (e.g., Tanaka et al. 2008; see also Braglia et al. 
2009 for the most massive clusters at z ~ 0.3). This could be the 
case of XMMU J0338. 7+0030, according to the projected distri- 
bution of its candidate member galaxies (Figs. 4 and 16). 

The spectroscopically confirmed, optically-selected proto- 
cluster CI 0332-2742 at z = 1.6 provides evidence of the build- 
up of the red sequence during its assembly (Kurk et al. 2009). 
Conversely, the X-ray luminous cluster XMMU J2235. 3-2557 
at z = 1.393 (Strazzullo et al. 2010) exhibits a tight red- 
sequence of massive galaxies, with overall old stellar popula- 
tions, generally early-type morphology, typically showing early- 
type spectral features and rest-frame far-UV emission consistent 
with very low SFRs. The massive cluster XMMU J1230.3+1339 
at z — 0.975 exhibits a richly populated red sequence, but 
also a few, spectro-photometrically classified, early-type galax- 
ies that potentially host some residual star-formation activity 
(Lerchster et al. 2011). Consistently, the X-ray selected clus- 
ter RDCS J09 10+5422, exhibits a red sequence although its 
early-type galaxy population appears to be still forming (Mei 
et al. 2006a; Tanaka et al. 2008). With a total mass of about 
10 14 M (Tanaka et al. 2008), RDCS J09 10+5422 is as mas- 
sive as XMMU J0338. 7+0030 and four times less massive than 
XMMUJ1230.3+1339 at z = 0.975 (Fassbender et al. 2011a; 
Lerchster et al. 201 1). 

If XMMU J0338.7+0030 lacks passively evolving galaxies, 
it hosts a high percentage of star-forming systems (late-type 



13 The observed red-sequences and colour thresholds for i' - Ks se- 
lected EROs (Pozzetti & Mannucci 2000 1 and J - Ks selected distant 
red galaxies (DRGs, Franx et al. 2003) reproduced in the two figures 
were corrected for the different filter transmission functions, assuming 
a 4 Gyr-old SSP with solar metallicity (see Pierini et al. 2004a). 



galaxies or starbursts), judging from the spectro-photometric 
classification of its candidate members (Fig. 15b). Nine of its 
10 candidate members that are spectro-photometrically classi- 
fied as starbursts are located well beyond a cluster-centric dis- 
tance of 1' (right panel of Fig. 4). Consistently, the spectro- 
scopically confirmed member at a cluster-centric distance of 2.3' 
(ID=12) is a bright galaxy with a quite strong [Oil] line emis- 
sion (Sect. 2.4.2). The other bright spectroscopic member with 
evidence of [Oil] line emission (ID=15) is a photometric mem- 
ber as well but is only 50" away from the cluster centre and sits 
on the red sequence in Fig. 15 as well as the two spectroscopic 
members classified as elliptical galaxies (ID=2 and ID=15) at 
cluster-centric distances of 28" and 2.1", respectively. 

Interestingly, there is evidence of sustained star-formation 
activity at the bright end of the galaxy luminosity function in 
the X-ray luminous cluster XMMU J 1007 .4+ 1237 at z = 1.555 
(Fassbender et al. 2011b) and in dense environments at z ~ 1 
(e.g., Elbaz et al. 2007; Gerke et al. 2007; Cooper et al. 2008). 
Several [O II] line emitters are found in the luminous, red galaxy 
population of XMMU J0302.2-0001 at z = 1.185 with a total 
mass of a few times 10 14 Mq (Suhadaet al. 201 1). Spectroscopic 
cluster members with [Oil] line emission can populate the lo- 
cus of old, passively evolving galaxies or be even redder than 
the red sequence (see XMMU J2235. 3-2557; Lidman et al. 
2008). Consistently, a non negligible fraction (20%) of the lu- 
minous infrared galaxies in 16 optically-selected galaxy clusters 
at 0.4 < z < 0.8 sits on the corresponding red sequences (Finn 
et al. 2010), the remaining 80% laying below. This is no sur- 
prise when the effect of dust attenuation on the observed SED 
of a galaxy is computed in a physical way, tailored to its star- 
formation mod43(cf. Pierini et al. 2004a, 2004b, 2005). 

In addition, star-forming galaxies at z ~ 1-2.5, selected in 
the rest-frame ultraviolet (UV) and with stellar masses of 10 10 - 
10 11 M0, exhibit evidence of complex configurations of stars 
and dusty interstellar medium (e.g., Noll et al. 2007, 2009). 
This points to the existence of dusty gas outflows, likely driven 
by feedback from supernovae, at high redshifts (cf. Steidel et 
al. 2010), but gravitational interactions cannot be excludecQ 
Extended distributions of dust not associated with the main bod- 
ies of galaxies have been discovered in the nearby Stephan's 
Quintet (Natale et al. 2010) and M81 triplet of galaxies (Walter 
et al. 201 1). Walter et al. (201 1) highlight the importance of tidal 
stripping for the metal enrichment of the intergalactic medium 
at high redshifts in addition to outflows. Here we note that the 
additional cooling of the ICM on dust grains impacts on the X- 
ray scaling properties of groups and clusters of galaxies (e.g., da 
Silva et al. 2009; Natale et al. 2010). Furthermore, early preheat- 
ing of the ICM is likely responsible for the observed evolution 
of the X-ray scaling relations of galaxy clusters out to z ~ 1.5 
(Reichert et al. 2011). 

In conclusion, XMMU J0338. 7+0030 appears as a bound, 
low-mass system of galaxies at z = 1.1 that is caught in an early 
assembly phase and where the shaping of the bright end of the 
galaxy luminosity function and nurturing effects on galaxy evo- 
lution (see Boselli & Gavazzi 2006 for a review) can be sought. 



14 Assumptions on the relative distribution of dust and stars as well as 
on the extinction law seem to impact the estimates of SFR and stellar 
mass in a systematic way (e.g., Kiipcii Yolda§ et al. 2007). 

15 Hence, some complexity in the relation between bolometric in- 
frared luminosity associated with re-emission by dust at mid-IR through 
sub-millimeter wavelengths and colour temperature of the warm and 
cold dust components can be expected. This could indeed have already 
been observed (cf. Hwang et al. 2010). 
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A full understanding of the properties of the galaxy population 
and ICM of this cluster requires a more in-depth study based on 
spectroscopy and multi-wavelength imaging, which offer better 
statistics and/or superior quality. 

5. Conclusions 

This paper presents, in particular, results from the first pointed 
observations of a high-z cluster with GROND, the seven-channel 
imager mounted at the MPI/ESO 2.2m telescope of La Silla, 
Chile. The target was the X-ray selected, weak, marginally 
extended source XMMU J0338.7+0030 in the XMM-Newton 
Distant Cluster Project (XDCP) survey. 

Follow-up imaging with OMEGA2000 at the 3.5m CAHA 
telescope enabled a comparison of the reddest z - H colours of 
the galaxies within a cluster-centric distance of 45" with predic- 
tions from simple stellar population evolutionary models, which 
suggested a redshift of z = 1.45 + 0.15 (lo - ). Later available 
VLT/FORS2 spectroscopy identified XMMU J0338.7+0030 as 
a galaxy cluster at z = 1.097 + 0.002 (lo - ) with four spec- 
troscopic members. From the determination of the flux in the 
soft 0.5-2 keV band, the bolometric luminosity of the cluster 
was inferred. The ensuing X-ray estimate of the total mass of 
XMMU J0338. 7+0030, based on a luminosity scaling relation, 
is equal to M 2 oo ~ 10 14 M©. 

About six hours of simultaneous imaging in the GROND 
g',r',i',z', J,H, Ks bands enabled us to extract photometry for 
832 sources detected down to z' AB ~ 26 (lcr). Existing coverage 
of the same region by the shallower Sloan Digital Sky Survey 
and 2 Micron All Sky Survey allows a robust photometric cali- 
bration to be obtained for the optical and near-infrared (IR) chan- 
nels of GROND, respectively. This is the basis for the applica- 
tion of the photometric redshift technique, as implemented in the 
publicly available code le Phare. 

With the GROND data, we confirmed the existence of an 
optical counterpart to XMMU J0338. 7+0030 on the basis of 
five consistent and independent pieces of evidence. In partic- 
ular, the Ks-band galaxy number counts for the cluster region 
imaged with GROND exhibit an excess by a factor of two 
at 20.5-22.5 AB mag with respect to those from survey ar- 
eas/deep fields. In addition, the photometric redshifts of the 
three spectroscopic members imaged with GROND (1.12 + 0.09, 
lcr) is consistent with the spectroscopic redshift of the clus- 
ter. Furthermore, the distribution of the 44 candidate members 
of XMMU J0338.7+0030 (with photometric redshifts within the 
range 1.01-1.23) in the i' - z' vs z' colour-magnitude diagram 
exhibits a red locus and a bluer one that are consistent with the 
red sequence and the distribution of the star-forming galaxy pop- 
ulation of the similarly massive cluster RDCS J09 10+5422 at 
Z= 1.106. 

XMMUJ0338.7+0030 seems to host a galaxy population 
that can still undergo significant bursts of star-formation activity, 
as also found in other XDCP high-z clusters. Part of these star- 
forming galaxies can exhibit colours that are at least as red as 
those of passively evolving galaxies, owing to dust attenuation, 
in agreement with results in the literature. A comparison of the 
projected distributions of the X-ray emitting plasma and the pho- 
tometric cluster members suggests that XMMU J0338. 7+0030 is 
likely caught in an early phase of its assembly: its star-forming 
galaxies fall in and face an environment that rapidly changes. 

Finally, we acknowledge that the identification of three spec- 
troscopic members out of the 832 GROND sources was essen- 
tial to determine the photo-z of the cluster as well as to select 
its candidate (photometric) members. However, it is remarkable 



that these two estimates of the redshift of XMMU J0338. 7+0030 
coincide within an uncertainty of 0.09, given that the so-called 
4000 A break (a fundamental feature for photo-z determination) 
falls in between the GROND i' and z' channels for a z = 1.1 
galaxy and the present unoptimized observations with GROND 
yielded a limited photometric accuracy in the near-//? bands. We 
also stress that no training of the solutions of the photo-z code 
was performed, because the available spectroscopic information 
is limited to only 0.6% of the GROND sources. 

For these reasons, this unique imager is useful for establish- 
ing optical counterparts to X-ray selected clusters, even in ab- 
sence of spectroscopic priors. This is particularly true for high 
throughput surveys for clusters like the eROSITA all-sky survey, 
which is expected to mostly deliver massive clusters at interme- 
diate redshifts. These clusters should host a wealthy population 
of old, passively evolving galaxies, which define a neat red se- 
quence in at least one of the colour-magnitude diagrams that can 
be built out of the GROND optical/near-IR photometry obtained 
in less than two hours of observations. 
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Fig.A.l. Distribution of the 68% range of the photo-z solu- 
tions (Az) as a function of the associated best-fit photo-z for 
the 436 GROND sources identified as galaxies (empty circles). 
Candidate members of XMMU J0338.7+0030 are reproduced 
with red (green) filled circles if they are within (outside) the bona 
fide cluster region; the three spectroscopic members among them 
are marked with blue filled circles (see Sect. 3.3.2). The dotted, 
short-dashed, long-dashed and solid lines represent values of Az 
equal to 0.4, 0.2, 0.1 and 0.05 x (1 + z), respectively. 
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Appendix A: Photometric redshifts: accuracy and 
reliability 

Here we discuss the accuracy of the photometric estimates of 
the distances to the individual sources in the 3.9 x 4.3 arcmin 2 
region of XMMU J0338.7+0030 imaged with GROND that is 
achieved through the GROND multi-band photometry and the 
photo-z code le Phare (see Sect. 3.2). At variance with works on 
survey areas/deep fields, we do so on the basis of photometric 
redshifts, mostly, because spectroscopic ones are available for 
only 0.6% of the GROND photometric sample (Sect. 2.4.1). 

In particular, we reproduce the behaviour of the 68% range 
of the photo-z solutions yielded by le Phare, Ai^l as a func- 
tion of the best-fit photo-z or z'-band magnitude for the 436 
sources out of 832 that are classified as galaxies by the best-fit 
solution (Figs. A.l and A. 2, respectively). The CFHTLS galaxy 
templates without emission lines are used. The 42 galaxies se- 
lected as fiducial cluster members (i.e., with 1.01 < z < 1.23, 
see Sect. 3.3.1) are colour-coded as in Fig. 15 (Sect. 3.2.2). 

Overall, the distribution of these photometric redshifts is 
well described by a straight line of equation Az ~ 0.05 x (1 + z) 
uptoz ~ 1.5; only 2 photo-z estimates (i.e., 0.4%) behave as out- 
liers, in the sense that they exhibit a value of Az > 0.4 x (1 + z). 
Unsurprisingly, photometric redshifts lower than 1 are assigned 

16 Given its definition, Az corresponds to twice the value expected for 
a Gaussian photo-z distribution centred at the best-fit value. 



Fig.A.2. Distribution of Az as a function of the z'-band mag- 
nitude for the same GROND sources displayed in Fig. A.l. 
Symbols have the same meaning as well. 




0.5 1 

Z SD3S 



Fig. A.3. Comparison of photometric redshifts for the 24 galax- 
ies that are in common between our GROND photometric cata- 
logue and the SDSS DR7 one (Sect. 2.4.1). For the SDSS photo- 
z's, uncertainties correspond to +lcr values. The red line repre- 
sents equality. 



with an accuracy that is better than Az ~ 0.05 x (1 + z), whereas 
photometric redshifts between 1 .5 and 3 exhibit a mean accuracy 
equal to Az ~ 0. 1 x ( 1 + z) but with a large scatter. In addition, an 
accuracy of 0.2 or better characterizes the photometric redshifts 
attributed to galaxies that are brighter than z' ~ 23.5, which cor- 
responds to a lOcr detection (see Sect. 2.3.2). 

The candidate member galaxies of XMMU J0338.7+0030 
tend to exhibit the lowest values of Az among the galaxies with 
the same z'-band magnitude (Fig. A. 2). However, three fiducial 
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Fig. A.4. Comparison of photometric and spectroscopic redshifts 
for the five sources in the GROND photometric catalogue ob- 
served with FORS2 (Sect. 2.4.1). For each spectroscopic value, 
we show the 68% range of the photometric estimates. The red 
line represents equality. 



Figure B.l shows the z' - H vs H colour-magnitude di- 
agram of the same cluster obtained from the GROND simul- 
taneous multi-band imaging, where total magnitudes are com- 
puted following method B (see Sect. 2.3.2). Its comparison with 
the analogous diagram in Fig. 3 (Sect. 2.2.2) reveals that the 
spectroscopic members with ID=2 and ID=16 exhibit z' - H 
colours that are about 0.8 Vega mag and 0.6 Vega mag bluer than 
their z - H colours listed in Table 1 (Sect. 2.2.2). Therefore, 
they lie closer to the spectroscopic member with ID=15 (with 
z' - H — 2.24 Vega mag, in full agreement with its z - H colour 
listed in Table 1) in Fig. B.l than in Fig. 3. In addition, the ex- 
istence of a red sequence is suggested in Fig. B.l. If we com- 
pute total magnitudes following method A, the spectroscopic 
members with ID =2 and ID=16 exhibit z' - H colours that are 
about 0.4 Vega mag bluer than their z - H colours in both cases, 
whereas the spectroscopic member with ID=15 exhibits a z' - H 
colour that is about 0.3 Vega mag redder than its z-H colour this 
time. The H-band magnitudes of these three galaxies can differ 
by up to +0.5 Vega mag with respect to those listed in Table 1, 
as a function of object and method (A or B). 



members within a cluster-centric distance of 1' exhibit values of 
Az between 0.3 and 1.1: one is an outlier (Fig. A.l). 

As a consistency check of our set of photometric redshifts, 
we compared photo-z's for 24 GROND sources with coun- 
terparts in the SDSS DR7 (see Sect. 2.3.2). The comparison 
SDSS photometric redshifts were calculated using a Neural 
Network method (Oyaizu et al. 2008). The agreement between 
these two independent sets of photo-z estimates is good up to 
z ~ 1. As expected, our values have significantly lower uncer- 
tainties since they are based on medium-deep photometry in 
the g',r',i',z', J, H, Ks bands (e.g., Bolzonella, Miralles, Pello 
2000). However, the u-band photometry that is available from 
the SDSS shallower imaging in the u, g,r, i, z bands (Abazajian 
et al. 2009 and references therein) better helps constraining 
photo-z's of galaxies at z < 0.4 (Bolzonella, Miralles, Pello 
2000). 

Finally, we reproduce photometric redshifts vs spectroscopic 
ones for the five extragalactic sources in the GROND pho- 
tometric catalogue that were observed with VLT-FORS2 (see 
Sect. 2.4.1) in Fig. A.4. In spite of a slight tendency towards 
higher values, our photometric redshifts, computed without any 
training of the solutions of the photo-z code le Phare, agree with 
the spectroscopic ones within +0. 1 . 



The method adopted to compute total magnitudes certainly 
impacts the estimated brightness of sources like these three 
galaxies, which are at the limit bewteen being considered partly 
resolved or point-like for resolutions of 1.1-1.3". However, it is 
likely the stacking of individual frames without PSF-matching in 
the case of the OMEGA2000 images that explains the previous 
differences in colours. 



Appendix B: On the z - H vs H colour-magnitude 
diagram of XMMU J0338.7+0030 

Here we investigate the origin of the discrepancy between the 
spectroscopic redshift of XMMU J0338.7+0030 (z = 1.097 ± 
0.002, Sect. 2.4.2) and the value that was suggested from the 
analysis of the z - H vs H colour-magnitude diagram obtained 
from the CAHA imaging (z ~ 1.45 ± 0.15, Sect. 2.2.2). There 
total magnitudes correspond to Kron aperture magnitudes i.e., 
computed as SExtractor MAG_AUTO (Sect. 2.1.2). 



We note that an SSP model with Salpeter IMF and solar 
metallicity exhibits z - H ~ 2.4 Vega mag when observed at 
z = 1.1, for a formation redhift between 3 and 10 (fig. 3 in 
Fassbender et al., 2011c). This colour is about 0.8 Vega mag 
bluer than the z - H colour of the same model observed at 
Z = 1.45 (Fig. 3). It also agrees well with the distribution of the 
spectro-photometrically classified elliptical galaxies in Fig. B.l, 
modulo (small) transformations between different z,H filters. 
At the same time, it is not inconsistent with the distribution of 
the spectroscopic members with absent/weak [Oil] line emission 
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Fig. B.l. z' - H vs H colour-magnitude diagram for all GROND sources in the XMMU J0338. 7+0030 region with flux detections 
in the z',H bands at a significance level > lcr, whatever the classification of the source by Le Phare (i.e., stars are included). 
Photometric errors (lcr) are also shown. As in Fig. 15, candidate members of XMMU J0338. 7+0030 are reproduced with red 
(green) filled circles if they are within (outside) the bona fide cluster region (panel a). The three spectroscopic cluster members in the 
GROND photometric catalogue are marked with blue filled circles there. Green asterisks represent the two QSO at a cluster-centric 
distance greater than V . Alternatively (panel b), photometric cluster members classified as galaxies are colour-coded according to 
their spectro-photometric type, which is E (red), Sbc(yellow), Scd (green), Im (cyan) or SB (blue). In each panel, the red short- 
dashed, solid and long-dashed lines represent the lcr, 3<x, and 5<x flux thresholds in the z',H bands (cf. Sect. 2.3.2). Photometric 
members classified as elliptical galaxies define a narrow red locus in this colour-magnitude diagram. 

(ID=2 and ID=15) in Fig. 3, taking into account the photometric 
uncertainties. 
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